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[Article  by  LM,  Arshavskiy,  C.Tech.S.,  V.A.  Kazakov  and  S.B.  Korolev,  engineers,  and  A.Ye.  Kroshilin, 
D.Tech.S.,  Nuclear  Power  Plant  SRI  [A-U]] 

[Text]  Analysis  of  the  dynamic  regimes  that  occur  in  the  process  of  operating  a  power  generating  unit 
shows  that  some  time-dependent  processes  involving  process  equipment  failures  can  lead  to  serious 
perturbations  in  certain  parameters  with  subsequent  generating  of  a  shutdown  signal.  The  turbogenerator 
load  shedding  regime  is  one  such  process,  wherein  one  of  the  two  turbine  feedwater  pumps  and  two  of  the 
four  main  circulating  pumps  which  operate  under  nominal  power  levels  are  disconnected.  The  primary 
reason  for  a  shutdown  occuring  in  these  regimes  is  that  reactor  load  reduction  is  too  slow.  Using 
accelerated  load  reduction  for  the  power  unit,  which  involves  dropping  a  cluster  of  control  rods  into  the 
core,  allows  us  to  avoid  a  shutdown,  thereby  improving  the  power  plant’s  economic  indicators,  reducing 
the  amplitude  of  thermal  interactions  with  the  equipment,  and  reducing  the  probability  of  a  dangerous 
radiation  leak  situation  as  the  reactor  makes  excursion  to  a  critical  state. 

We  present  the  results  of  a  numerical  and  theoretical  analysis  of  the  efficacy  of  an  accelerated  load 
reduction  system  (ALR).  We  used  the  KIPR  integrated  mathematical  model  ^  as  a  tool  in  our 
investigation.  We  look  at  three  separate  classes  of  nuclear  power  plant  operating  regimes  in  which  an 
accelerated  load  reduction  system  is  employed. 

Dynamic  regimes  during  turbogenerator  load  shedding.  In  regimes  where  turbogenerator  load  is  shed 
(down  to  no-load,  internal-needs  only,  or  zero)  a  strong  increase  in  primary  feedwater  inlet  pressure  is 
observed,  which  leads  to  full  opening  of  all  the  BRU-A  atmospheric  steam  dump  bypass  valves. 

According  to  the  design  algorithm,  in  reactor  load  reduction  one  may  note  the  following  negative 
phenomena  that  accompany  this  class  of  transient: 

•  if  certain  of  the  steam  dumps  fail  (either  the  BRU-K  into  the  condenser  or  the  BRU-A  into  the 
atmosphere)  the  reactor  may  be  shut  down  because  of  the  pressure  in  the  steam  generator; 

•  if  the  BRU-A  is  not  closed,  the  reactor  installation  may  be  cooled  at  an  unsafe  rate; 

•  an  increase  in  pressure  in  the  secondary  loop  causes  an  increase  in  pressure  in  the  primary  loop, 
resulting  in  opening  of  all  the  spray  nozzles,  which  is  the  cause  of  a  significant  drop  in  pressure  in  the 
primary  loop  during  the  final  phase  of  a  transient.  If  the  spray  nozzles  are  not  activated,  a  scram  signal 
is  generated; 
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•  a  strong  rise  in  temperature  on  the  cold  side; 

•  a  high  rate  of  steam  flow  into  the  condenser  via  the  BRU-K  steam  dump  may  cause  the  turbine 
feedwater  pump  to  be  shut  down. 

Employing  ALR  in  these  regimes  may  either  eliminate  some  of  these  phenomena,  or  at  least  significantly 
lessen  their  negative  effect  on  the  operation  of  the  power  unit.  This  is  confirmed  by  the  results  of  several 
numerical  studies  of  flows  in  these  classes  of  regime  (taking  into  account  possible  failures  during  operation 
of  the  equipment)  comparing  the  use  of  regular  algorithms  and  those  using  ALR. 

We  present  the  results  of  some  of  these  analyses  (see  Fig.  1).  One  regime  corresponds  to  the  load 
dropping  to  zero  upon  failure  of  two  of  the  BRU-K  steam  dumps.  The  initial  state  of  the  unit 
corresponds  to  the  end  of  a  refueling  cycle  with  a  power  level  of  102%.  This  is  confirmed  by  the  Figure. 

We  will  look  now  at  the  dynamics  of  a  transient  process  as  turbogenerator  load  is  shed  down  to  a  no-load 
or  idle  condition  using  ALR  at  a  nominal  power  level,  which  was  realized  on  7  Jun  87  on  the  No  3  unit  of 
the  Zaporozhe  AES.  The  experimental  and  numerical  results  of  this  dynamic  test  are  shown  in  Fig.  2. 

By  plunging  Cluster  I  into  the  core,  in  10  s  the  neutron  power  was  52%  of  (54%  of  A'nom 

calculation).  Further  reactor  load  reduction  down  to  40%  was  carried  out  by  a  power  limiter  acting  on 
Control  Rod  Cluster  X. 


A  graph  showing  variation  in  pressure  in  the  primary  feedwater  inlet  is  given  in  Fig.  2.  The  BRU-K 
steam  dump  valves  began  opening  in  proportion  to  the  amount  of  load  shed  and  were  fully  open  after 
18  s.  The  valves  started  to  close  at  40  s  (in  the  experiment,  44  s)  in  step  with  the  pressure  drop  in  the 
secondary  loop.  In  the  experiment  one  of  the  BRU-K  steam  dumps  failed  to  open  as  it  should,  and  this 
was  allowed  for  in  the  calculations.  The  maximum  pressure  in  the  in  the  primary  feedwater  inlet  was  7.0 
MPa,  and  7.1  MPa  in  the  steam  generators,  which  is  lower  than  that  needed  to  trigger  atmospheric 

exhaust  via  the  BRU-A.  In  the  experiment  one  of 


the  BRU-As  did  open  slightly,  by  5 — 10%.  We 
should  note  that  the  actual  flow  characteristics  of 
the  steam  dumps  may  differ  strongly  from  design. 
Thus  an  analysis  of  this  regime  shows  that  the 
pressure  in  the  feedwater  inlet  stabilized  at  a  level  of 
roughly  5.9  MPa  at  a  power  level  of  39%  of  A'nom 
and  a  feedwater  flowrate  of  about  40%  of  G^nom’ 
with  three  of  the  BRU-K  valves  40%  open.  The 
feedwater  temperature  equalized  at  about  170°  C. 
Simple  calculations  show  that  the  steam  flowrate  in 
such  conditions  through  a  fully-open  BRU-K  is 


1.20 — 1.35  of  the  design  values  (the  lack  of  precision 
in  determining  the  flowrate  is  explained  by 
incompleteness  of  data  about  flowrates  through  the 
BRU-SN  steam  dump  and  on  to  the  turbine  as  load 
is  shed  to  an  idle  or  no-load  condition,  as  well  as 

Figure  L  Graph  of  Variation  in  Parameters  in  the 
Closed  SRK  regime  (End  of  Refueling  Cycle,  Minus 
Two  BRU-K  Steam  Dumps):  i -Neutron  power;  2  and 
5-pressure  in  main  feedwater  inlet  and  pressurizer;  4~ 
temperature  on  the  cold  side  (solid  line,  ALR  in  use;  dashed 


line,  ALR  not  in  use) 


2 


This  analysis  shows  that  using  ALR  and  opening 
the  BRU-K  in  proportion  to  the  load  shed  in  this 
regime,  even  when  one  of  the  BRU-K  fails,  will 
ensure  that  enough  pressure  is  maintained  to  keep 
the  BRU-A  closed  (one  BRU-A  opened  by  a  total 
of  only  5—10%)  and  allows  us  to  avoid  opening  up 
the  pressurizer  spray  nozzles  during  the  transient. 

The  good  agreement  between  experimental  and 
calculated  results  confirms  that  the  KIPR 

Figure  3.  Calculated  Curves  of  Variation  in 
Pressure  in  Pressurizer  (a)  and  Main  Feedwater 
Inlet  {h)  in  the  Closed  SRK  Regime  Using  ALR:  i- 

Start  of  a  refueling  cycle,  all  steam  dumps  operating;  S-End 
of  a  refueling  cycle,  all  steam  dumps  operating;  5-Start  of  a 
refueling  cycle,  minus  four  of  the  BRU-K  steam  dumps;  4- 
end  of  a  refueling  cycle,  minus  four  BRU-K  steam  dumps 


Figure  1.  Graph  of  Variation  in  Parameters  in  the 
Load  Shedding  Regime  Down  to  No-Load  at  the 
Zaporozhe  AES:  i -Neutron  power;  2  and  5-pressure  in 
main  feedwater  inlet  and  pressurizer;  ^-temperature  on  the 
cold  and  hot  sides;  ^-mean  temperature  in  the  primary  loop 
(solid  line  is  experiment,  dashed  line  is  calculation) 


incompleteness  of  data  on  the  degree  of  nonlinearity 
of  the  BRU-K  valve  flowrate  characteristic  curve). 
When  performing  the  calculations  for  this  regime  the 
flowrate  through  the  open  BRU-K  was  taken  to  be 
1-25  Gdesign  =  1125  t/h. 
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integrated  software  package  is  accurate  to  a  satisfactory  degree,  and  lets  us  carry  out  numerical  analysis 
of  a  number  of  other  possible  variants  in  the  flows  of  other  regimes  that  employ  ALR. 

The  following  variants  were  examined: 

1.  Load  shedding  down  to  zero;  start  of  a  refuel  cycle;  all  steam  dumps  in  operation. 

2.  Load  shedding  down  to  zero;  start  of  a  refuel  cycle;  four  of  the  BRU-K  steam  dumps  fail. 

3.  Load  shedding  down  to  zero;  end  of  a  refuel  cycle;  all  steam  dumps  in  operation. 

4.  Load  shedding  down  to  zero;  end  of  a  refuel  cycle;  four  of  the  BRU-K  steam  dumps  fail. 

In  carrying  out  the  calculations  the  flowrate  through  the  steam  dumps  was  taken  to  be  1125  t/h  at  a 
pressure  of  6.28  MPa.  A  curve  of  the  variation  in  pressure  in  the  primary  and  secondary  loops  in  this 
regime  is  shown  in  Fig.  3.  From  analysis  of  the  curve  we  may  draw  the  following  conclusions: 

•  at  the  beginning  of  a  refuel  cycle  no  opening  of  the  spray  nozzles  in  the  pressurizer  was  observed.  Thus, 
even  when  four  of  the  BRU-K  steam  dumps  fail  the  maximum  pressure  in  the  primary  loop  exceeds  the 
initial  pressure  by  only  0.28  MPa,  and  does  not  reach  a  level  where  the  spray  nozzles  are  operated.  At  the 
end  of  a  cycle  (in  variant  4)  both  spray  nozzles  are  opened  (by  84  and  50%  respectively)  and  the 
maximum  pressure  in  the  pressurizer  is  16.35  MPa  at  20  s.  However,  if  they  fail  this  does  not  lead  to  a 
shutdown;  the  pressure  must  increase  to  17.3  MPa  for  this  to  happen. 

•  in  all  regimes  the  minimum  pressure  in  the  primary  loop  is  no  lower  than  14.5  MPa; 

•  at  the  beginning  and  at  the  end  of  a  refuel  cycle,  with  all  steam  dumps  operating,  the  pressure  in  the 
main  feedwater  inlets  is  not  more  than  1.25  MPa.  There  was  a  momentary  opening  of  the  BRU-A 
observed;  in  variant  1,  by  26%  and  in  variant  3,  by  34%.  We  note  that  if  the  initial  pressure  in  the 
feedwater  inlet  was  5.88  MPa  (60  kgf/cm^),  which  is  what  occured  in  the  dynamic  regime  we  looked  at 
earlier,  then  there  is  no  opening  of  the  BRU-A; 

•  when  all  four  BRU-K  fail,  all  the  BRU-A  open  fully.  The  maximum  pressure  in  the  main  feedwater 
inlet  (and  in  the  steam  generators,  since  in  these  variants  they  are  approximately  the  same)  is:  in  regime 
2,  2.68  MPa  and  in  regime  4,  7.74  MPa,  which  is  quite  near  shutdown  pressure  level  for  the  steam 
generator  of  7.85  MPa. 

Dynamic  regimes  with  one  of  two  operating  turbine  feedwater  pumps  disconnected.  Disconnection  of  one 
of  the  two  operating  turbine  feedwater  pumps  results  in  serious  perturbations  in  the  feedwater  flow  to  the 
steam  generator,  and  thus  the  primary  problem  in  this  process  is  to  maintain  a  nominal  water  level  in  the 
steam  generator. 


At  nominal  power  levels,  a  shutdown  results  if  this  regime  runs  its  course  (either  from  the  500  mm  drop 
in  level  in  one  of  the  steam  generators,  or  from  the  fact  that  two  of  the  operating  main  circulation  pumps 
are  disconnected,  causing  a  200  mm  drop  in  level  in  the  corresponding  steam  generators).  Accelerated 
load  reduction  of  the  unit  makes  it  possible  to  rapidly  reduce  the  thermal  flux  to  the  steam-water  mixture 
in  the  steam  generator  from  the  primary  loop  side,  and  thus  to  lower  the  steam  flowrate  from  the  steam 
generator  and  reduce  the  imbalance  in  feedwater  and  steam  flowrates,  which  is  the  primary  reason  for  the 
sharp  drop  in  water  levels  in  the  steam  generators. 

However,  dropping  a  rod  cluster  in  this  regime  results  in  a  pressure  drop  in  the  main  feedwater  inlet,  if 
the  EGSR  turbine  regulator  is  in  the  operating  regime  for  maintaining  a  nominal  steam  pressure  before 
the  turbogenerator  SRK  valve  throughout  the  transient  process.  Thus  we  need  a  different  algorithm  for 
reducing  the  load  on  the  turbine.  Initially  we  proposed  the  following: 
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•  reducing  the  load  on  the  turbine  by  150  MW  in  the  first  1  s;  and 

•  further  reducing  the  load  down  to  500  MW  at  a  rate  of  15  20  MW/s. 

Analysis  of  the  regime  in  which  one  of  two  operating  turbine  feedwater  pumps  is  disconnected,  and  ALR 
is  employed,  is  carried  out  using  the  experimental  and  numerical  results  from  two  dynamic  tests 
performed  on  VVER- 1000  power  generating  units. 

Regime  1  was  carried  out  on  28  Jan  89  on  the  No  2  unit  of  the  South  Ukraine  AES  at  the  end  of  a 
refueling  cycle  (270  eff.  days). 

Regime  2  was  carried  out  on  29  May  87  on  the  No  3  unit  of  the  Zaporozhe  AES  at  the  start  of  a  refueling 
cycle  (127  eff.  days). 

The  initial  algorithm  for  turbine  load  reduction  was  employed,  and  in  the  experiment  the  EGSR  regulator 
was  in  the  “process  safety”  mode  when  after  1  s  the  turbine  load  was  reduced  by  160  MW  and  thereafter 
at  a  rate  of  10  MW/s  down  to  a  level  of  500  MW  (38  s  into  the  process),  after  which  the  turbine  went 
into  the  RD-1  regime  with  a  regulator  setting  of  6.0  MPa.  The  inital  turbogenerator  load  reduction  led  to 
an  increase  in  the  pressure  at  the  main  feedwater  inlet  and  the  activating  of  the  BRU-K  valves.  Three 
BRU-K  valves  were  20%  closed  15  s  into  the  transient.  The  pressure  in  the  feedwater  inlet  reached  its 
maximum  value  of  6.45  MPa  at  26  s,  after  which  it  fell  in  proportion  with  the  thermal  power  of  the  unit 
and  the  operation  of  the  steam  dumping  equipment. 

As  has  already  been  noted,  defining  the  EGSR  turbine  regulator  operating  algorithm  is  the  most 
important  step  in  implementing  ALR.  The  initial  algorithm  leads  to  an  increase  in  pressure  in  the  main 
feedwater  inlet  and  opening  of  the  BRU-K,  which  is  undesirable  because  of  possible  failures  of  these 
regulators.  A  sharp  drop  in  pressure  at  the  feedwater  inlet  is  also  inadmissible.  The  best  thing  would  be 
a  turbine  load  reduction  algorithm  such  that  the  pressure  at  the  feedwater  inlet  is  maintained  somewhere 
in  the  range  p^^^±0.2  MPa. 

An  optimal  turbogenerator  load  reduction  algorithm  in  a  regime  where  one  of  the  two  operating  turbine 
feedwater  pumps  is  shut  off  was  selected  on  the  basis  of  a  multivariate  calculation  using  the  KIPR 
software  package.  It  assumes  that  the  EGSR  regulator  operates  during  this  transient  in  a  regime  that 
maintains  steam  pressure  before  the  turbine  (RD-1),  however,  this  pressure  is  variable  and  is  determined 
from  the  graph  in  Fig.  4. 

This  algorithm  was  implemented  in  regime  2  (Fig.  5).  Within  1.4  s  after  disconnecting  turbine  feedwater 
pump  TFP-2  a  signal  was  sent  to  activate  power  limiting  and  drop  the  first  control  rod  cluster.  Turbine 
load  reduction  occurred  at  the  same  time  as  reactor  load  reduction.  Upon  a  signal  for  accelerated  load 
reduction  to  the  EGSR  the  RD-1  pressure  maintainance  regime  was  activated  and  within  4  s  that  steam 
pressure  began  to  change  in  accordance  with  the  selected  algorithm  (Fig.  4).  Turbine  load  reduction  was 
initiated  after  6  s  and  was  terminated  at  37  s  by  a  false  switching  of  the  automatic  output  power 
regulator  to  the  T  regime  and  of  the  EGSR  to  the  regime  where  turbine  power  is  kept  constant  at  660 
MW  (for  a  reactor  power  of  49%).  This  was  taken  into  account  in  performing  the  numerical  calculations 
for  this  regime.  The  imbalance  in  reactor  and  turbine  power  led  to  a  drop  in  pressure  at  the  main 
feedwater  inlet.  However,  the  good  agreement  between  the  numerical  and  experimental  pressure  behavior 
at  the  feedwater  inlet  up  to  the  37  s  point  allows  us  to  conclude  that,  were  it  not  for  the  false  switching  of 
the  EGSR  taking  it  from  the  RD-1  regime,  the  pressure  at  the  feedwater  inlet  could  have  been  maintained 
within  the  prescribed  limits.  Fig.  5  shows  the  numerical  behavior  of  the  pressure  at  the  feedwater  inlet 
for  precise  implementation  of  the  selected  algorithm. 
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Figure  4.  Graph  of  Variation  of  Pressure  Setting  in 
Turbine  EGSR  Regulator 

Figure  5.  Graph  of  Variation  in  Parameters  in  the 
Regime  Where  One  Turbine  Feedwater  Pump  is 
Disconnected  at  the  Zaporozhe  AES:  i -Electrical 
power;  2  and  5-Flowrates  of  feedwater  to  and  steam  from 
steam  generator;  ^-pressure  in  main  feedwater  inlet;  5-level 
in  steam  generator  (solid  line  is  experiment,  dashed  line  is 
calculated,  dash- dot  line  is  calculated  while  observing 
EGSR  algorithm) 


From  analysis  of  the  regime  where  one  of  two 
operating  turbine  feedwater  pumps  is  disconnected, 
we  may  conclude  the  following: 

•  ALR  in  this  regime  allows  us  to  avoid  shutdown 
because  of  reduced  water  levels  in  the  steam 
generator; 

•  the  selected  algorithm  for  EGSR  regulator 
operation  ensures  that  the  pressure  at  the  feedwater 
inlet  will  be  maintained  within  the  prescribed 
limits. 

Dynamic  regimes  when  two  of  the  four  operating 
main  circulating  pumps  are  disconnected.  Regimes 
in  which  various  numbers  of  the  main  circulating 
pumps  are  disconnected  involve  a  decrease  in 
coolant  flowrate  through  the  core.  The  absence  of 
DNB  over  the  fuel  rod  surface  is  the  criterion  for 
reliable  core  cooling.  The  output  power  limiter 
reduces  reactor  load  down  to  67%  to  meet 

this  criterion  when  one  of  the  four  main  circulating 
pumps  is  disconnected.  If  two  of  the  four  main 
circulating  pumps  are  disconnected  (either 
simultaneously  or  within  a  period  of  time  of  less 
than  70  s)  then,  as  designed,  a  shutdown  signal  is 
sent.  This  is  because  if  the  shutdown  signal  were 


not  sent,  DNB  could  arise  13  s  into  the  process  (according  to  calculations  using  the  Dinamika  program 
over  maximumally  thermal-stressed  fuel  rods  for  an  adverse,  albeit  possible,  combination  of  parameters. 
As  we  will  show  further  on,  using  accelerated  load  reduction  of  a  unit  in  the  regime  where  two  main 
circulating  pumps  are  powered  down  lets  us  avoid  DNB  and,  consequently,  the  need  for  a  shutdown. 


The  results  of  processing  the  data  from  dynamic  tests  where  two  of  four  operating  main  circulating  pumps 
are  disconnected  using  ALR,  as  carried  out  on  4  Jun  87.  at  the  No  3  unit  of  the  Zaporozhe  AES,  are  shown 
in  Fig.  6.  As  the  two  pumps  were  disconnected,  after  1.4  s  the  power  limiters  were  activated  and  a  signal 
was  sent  to  drop  the  first  cluster  of  scram  rods.  Twenty-five  seconds  into  the  transient  the  control  rods 
had  reduced  the  reactor  load  down  to  49%,  and  the  automatic  output  power  regulator  was  connected  to 
control  It.  As  a  result  of  the  increased  pressure  at  the  feedwater  inlet  (by  roughly  1.5  bar)  the  output 
power  regulator  was  switched  over  to  the  T  regime  (maintaining  the  pressure  in  the  secondary  loop)  and 


Figure  6.  Graph  of  Variation  in  Parameters  in  the 
Regime  Where  Two  Main  Circulating  Pumps  are 
Disconnected  at  the  Zaporozhe  AES:  i  and  2- 

Electrical  and  neutron  power;  5-Flowrate  through  the  core; 
4  and  5-pressure  in  the  pressurizer  and  the  main  feedwater 
inlet  (solid  line  is  experiment,  dashed  line  is  calculated) 

reactor  load  was  reduced  down  to  36%  of  A^nom' 
Upon  a  signal  from  the  ALR  to  the  EGSR  the  RD- 
1  regime  was  engaged  and  preset  variation  in 
pressure  as  per  Fig.  4  was  initiated.  This  EGSR 
operating  algorithm  ensured  that  the  pressure  in 
the  primary  and  secondary  loops  is  maintained 
within  permissible  limits. 

Disconnection  of  two  of  four  operating  main 
circulating  pumps  is  accompanied  by  a  reduction  in 
the  flowrate  through  the  core  (Fig.  6)  and  the 
danger  of  a  breakdown  in  reliable  core  cooling 
conditions.  However,  dropping  in  a  cluster  of  rods 
leads  to  a  sharp  decline  in  power  output  and  a 
reduction  in  the  thermal  flux  coming  off  the  fuel 
rod  surfaces.  As  a  result  7^,  the  safety  margin  to 
DNB,  defined  as  the  ratio  of  current  thermal  flux 
to  critical  thermal  flux  under  the  given  conditions 
for  a  maximally-thermal  stressed  fuel  rod,  remains 
greater  than  unity  over  the  course  of  the  entire 
process.  The  minimum  value  of  t]  is  1.49  at 


roughly  2.2  s  (the  stationary  value  was  1.525).  We  ^  ^ 

should  note  that  the  critical  thermal  flux  calculations  were  performed  according  to  the  methods  in  ’  ,  the 
shape  of  the  neutron  field  and  the  coefficient  of  nonuniformity  over  the  core  volume  were  taken  from  the 
calculations  of  the  Dinamika  software,  and  the  rest  of  the  parameters  are  experimental. 


In  order  to  obtain  a  conservative  estimate  of  the  minimum  safety  margin  to  DNB  t}  in  the  regime  where 
two  main  circulating  pumps  are  simultaneously  disconnected,  we  performed  calculations  for  the  most 
adverse  combination  of  unit  parameters.  The  initial  conditions  were  as  follows:  N  —  107%  of  A'nom’ 
^pressurizer  =  15.4  MPa  (157  kgf/cm^),  =  6.08  MPa  (62  kgf/cm^),  and  flowrate  through 

the  reactor  G  =  80,000  m^/h.  The  shape  of  the  neutron  field  and  the  coefficient  of  its  nonuniformity  were 
taken  from  and  the  delay  before  dropping  a  cluster  of  control  rods  was  1.65  s. 

The  results  of  the  appropriate  calculations  show  that,  for  this  kind  of  combination  of  parameters,  r]  for  a 
thermally-stressed  fuel  rod  remains  substantially  greater  than  unity.  The  minimum  value  of  the 
coefficient  was  1.36  at  2.1  s;  the  initial  value  was  1.394.  The  calculations  we  performed  have 
demonstrated  that  when  two  of  the  main  circulating  pumps  are  simultaneously  disconnected,  using  ALR 
ensures  reliable  reactor  core  cooling. 

We  examine  now  the  regimes  in  which  two  of  the  main  circulating  pumps  are  disconnected  one  after  the 
other.  When  the  first  is  disconnected  it  sets  off  the  PZ-1  signal,  and  the  power  limiters  begin  reducing 
reactor  load  by  dropping  a  control  rod  cluster  into  the  core.  After  the  second  main  circulating  pump  is 
disconnected  (if  the  interval  of  time  between  disconnection  of  both  pumps  is  less  than  70  s)  the  first 
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cluster  is  lowered  into  the  core  after  a  delay  of  1.65  s.  In  the  calculations  we  analyzed  the  regime  where 
two^  main  circulating  pumps  were  successively  disconnected  in  an  interval  of  between  5  and  70  s.  The 
initial  conditions  of  the  unit  and  the  neutron  physics  properties  of  the  core  correspond  to  that  described 
earlier.  From  the  standpoint  of  ensuring  reliable  core  cooling,  the  most  dangerous  regime  is  that  in  which 
the  delay  in  disconnecting  the  second  pump  is  At  =  18  s,  because  when  the  first  pump  has  been 
disconnected  and  the  power  limited  activated  the  minimum  value  of  the  coefficient  -q  is  reached  18  s  into 
the  transient.  However,  even  in  this  case  q  remains  substantially  greater  than  unity. 

In  the  At  =  10  s  regimes  no  excursion  of  the  unit  parameters  beyond  permissible  limits  was  noted,  with 
the  exception  of  a  combination  of  the  disconnected  pumps.  The  calculations  show  that  if  we  first 
disconnect  the  main  circulating  pump  for  the  loop  connected  to  the  pressurizer  (the  fourth  loop  in  an 
actual  unit),  then  the  shutdown  triggering  point  is  reached  on  account  of  the  drop  in  pressure  in  the 
primary  loop.  The  results  of  modelling  this  regime  are  given  in  Fig.  7.  The  principal  reason  for  the  drop 
in  pressure  in  the  primary  loop  to  the  shutdown  triggering  point  of  13.73  MPa  (140  kgf/cm^)  is  the 
massive  reactor  load  reduction  following  the  lowering  of  the  control  rods  as  a  result  of  a  PZ-1  signal  being 
generated  from  the  temperature  on  the  hot  side  (because  of  a  large  influx  of  saturated  water  from  the 
pressurizer  as  compared  to  underheated  water  from  the  VKS). 


This  adverse  effect  of  a  strong  increase  in  temperature  on  the  hot  side  is  manifest  only  for  this 
combination  of  connected  pumps  and  with  a  delay  in  disconnecting  the  main  circulating  pumps  of 
5  s  <  Ar  <  20  s.  For  Ar  <  5  s  the  flowrate  through  the  fourth  loop  at  the  moment  ALR  is  engaged  is 
adequate  and  the  inflow  of  coolant  from  the  pressurizer  does  not  lead  to  a  strong  increase  in  temperature 


on  the  hot  side.  Similarly,  for  Ar  >  20  s  a  reverse 
flow  of  coolant  in  these  loops  precludes  an  increase 
in  temperature  and,  consequently,  the  prophylactic 
safety  measures  that  are  the  cause  of  the  massive 
reactor  load  reduction. 

Calculations  for  the  Ar  =  60  s  regime  show  that 
ALR  is  engaged  the  moment  the  power  limiters 
reduce  the  reactor  load  down  to  a  power  level  of 
75—80%  of  nominal,  which  may  lead  to  setting  off  a 
shutdown  from  the  lowering  of  the  pressure  in  the 
primary  loop. 

A  numerical  and  experimental  analysis  of  the 
regimes  in  which  two  of  four  operating  main 
circulating  pumps  are  disconnected  using  ALR 
allows  us  to  draw  the  following  conclusions: 

•  an  operating  algorithm  for  the  EGSR  turbine 
regulator  selected  on  the  basis  of  these  calculations 


Figure  7.  Calculated  Curves  of  Variation  in 
Parameters  in  the  Regime  Where  Two  Main 
Circulating  Pumps  are  Disconnected  With  a  Delay 

of  10  s:  ; -Neutron  power;  ^-Flowrate  Through  DT  (G  >  0 
upon  leaving  pressurizer);  ^-pressure  in  pressurizer;  5- Hot 
side  temperature  (solid  line  is  actual,  dashed  line  is  taking 
into  account  delay  in  thermocouple  readings) 
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will  maintain  unit  parameters  within  permissible  limits,  which  is  confimrmed  by  the  results  of  dynamic 
tests  in  which  two  of  the  main  circulating  pumps  were  simultaneously  disconnected; 

•  reliable  core  cooling  is  ensured  in  all  regimes,  even  for  adverse  combinations  of  regulator  parameters; 

•  in  two  cases  two  of  the  main  circulating  pumps  are  disconnected  one  after  the  other  a  shutdown  is 
possible  due  to  an  increase  in  the  temperature  of  the  primary  loop,  which  causes  a  massive  reactor  load 
reduction  during  the  transient.  In  the  first  case,  in  which  the  delay  until  disconnecting  the  second  pump 
in  in  the  5 — 20  s  range  (and  the  pump  connected  to  the  pressurizer  is  disconnected  first),  this  is  associated 
with  the  temperature  on  the  hot  side  setting  off  a  PZ-1  signal.  In  the  second  case,  in  which  the  second 
pump  is  switched  off  within  50  to  70  s  after  disconnection  of  the  first,  this  is  associated  with  a  reduced 
neutron  power  and  pressure  in  the  primary  loop  at  the  moment  ALR  is  engaged.  In  all  other  cases, 
including  simultaneous  disconnection  of  pumps,  no  tripping  of  the  shutdown  was  noted. 

Thus  the  results  of  a  detailed  numerical  analysis  of  the  dynamics  of  a  power  generating  unit  in  regimes 
where  ALR  operates  shows  the  advisability  and  efficacy  of  its  use.  The  numerical  study  allows  us  to 
refine  unit  regulating  system  algorithms  and  to  study  the  dynamics  of  regulators  over  a  broad  range  of 
operating  regimes. 
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celikov  and  B.S.  Disik,  engineers  at  the  Ignalinsk  AES] 

uninterrupted  operation  of  the  process  water  supply  system  that  ensures  supply  of  water  to  the  orimarv 
and  secondary  entities  of  the  generating  unit,  including:  ^  ^ 

•  the  process  condensers; 

•  systems  for  cooling  bearing  assemblies  in  the  equipment; 

•  reactor  cooling  systems; 

•  coolant  loops  in  the  emergency  system  channels  of  the  reactor; 

•  heat  exchanger  apparatus  in  the  ventilation  and  water  treatment  systems,  etc. 

The  safe  operation  of  the  reactor  installation  in  an  AES  in  large  part  depends  on  the  reliable  operation  of 

all  subassemblies  in  the  process  water  supply  system  (PWSS),  in  particular  during  transient  and 
emergency  unit  operating  regimes^  ^  <.ransiem;  ana 

?!!  branched  hydraulic  network;  one  of  its  basic  features  is  the  differing  requirements 

imposed  by  the  different  consumers  in  terms  of  flowrates  and  water  pressure.  It  is  assumed  that  these 

Ziots  mC'in  'Hr  '  consumers  fluctuates  in  the  range 

irom  ZAi  to  10.28  m  /s.  In  addition,  there  are  seasonal  fluctuations  in  flowrates. 

i"  •  •’WSS  is  ths  pumping  equipmunt  that  supplies  .ate,  to  the 
H  Ensuring  the  reliability  .„d  proper  operation  of  a  pump  set  is  a  priority 

I"'*®  oft™  tlie  operating  conditions  are  examined  without  reference  to  its  connection 

tadi  mS::^”af  if  d  ’’"''Tf  of  -  individual  pumpt 

Hardly  monitored  at  all  during  the  monitoring  of  PWSS  parameters. 

^  1000V-y63  pump  set.  Analysis  of  the  results  of  our  observing  the  operation  of  pumps 
at  two  AES  power  units  shows  that  the  range  of  water  delivery  covers  the  entire  performance  range  of 

Sv“:n  s'uX  diitni"’  - 
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Figure  1.  Performance  Characteristic  Curve  for  a 
Centrifugal  Pump  in  a  PWSS 

Figure  2.  Graph  of  Vibration  in  Lower  Bearings  as 
a  Function  of  Operating  Regime:  i,  2,  3  are  2,  2.5  and 

4  Hz;  4  is  the  rotation  frequency;  5  is  twice  the  rotation 
frequency;  6  is  the  blade  frequency. 

If  we  look  at  a  pump  characteristic  curve  (see  Fig. 

1),  in  certain  cases  we  note  regime  points  that  go 
beyond  the  performance  limits.  This  becomes 
possible  because  the  useful  part  of  the  pump  head 
characteristic  curve,  especially  at  low  deliveries,  is 
gently  sloping  and  either  variations  in  the  operating 
regime  of  the  network  or  of  individual  differences  in 
pump  head  and  the  resistance  to  this  head  in  the 
lines  has  a  substantial  impact  on  the  location  of  the 
pump  regime  points  on  the  operational  characteristic 
curve.  In  turn  it  was  found  that  the  excursion  of  a 
pump  from  the  performance  range  is  characterized 
by  a  sharp  increase  in  dynamic  load  acting  on  the 
pump  set  housing  and  its  supports:  vibration 
increases,  as  does  the  intensity  of  hydrodynamic  pressure  pulsations.  It  was  found  that  pump  housing 
vibration  when  the  regime  departs  from  the  performance  range  may  be  triple  the  level  of  vibration 
relative  to  normal  operating  conditions. 

One  of  the  necessary  conditions  for  correct  operation  of  centrifugal  pump  sets  during  their  operation  in 
complicated  hydraulic  networks  is  the  timely  determination  of  the  quantitative  and  qualitative 
dependences  of  dynamic  loading  on  the  operating  regimes  of  the  pumps  themselves  and  the  PWSS. 

to  discover  these  dependences  we  studied  both  the  dynamic  characteristics  of  pumps  and  certain  of  the 
process  parameters:  flowrate  and  pressure  in  PWSS  pressurized  lines,  pressure  in  the  pump  discharges,  and 
the  current  draw  of  the  electric  drive  motors. 

Vibration  pickups  were  mounted  onto  the  intake  pipes,  the  lower  guide  bearings  and  the  drive  motor 
housing.  The  hydrodynamic  pressure  was  measured  in  the  region  below  the  pump  impeller,  in  the  transits 
of  the  spiral  case,  and  at  the  discharge  pipe.  The  measurements  were  performed  using  domestic  vibration 
measuring  equipment,  the  VI6-6TN  in  combination  with  a  DD-6S  pressure  pickup  and  DU-5S  vibration 
accelerometer,  as  well  as  some  equipment  from  Bruell  and  Kerr:  a  2120  spectrum  analyzer  and  2305 
printer.  A  study  was  made  of  three  operating  regimes  with  variation  of  flowrate  in  the  network  from  2.47 
to  3.17  m^/s.  The  operating  parameters  of  the  pump  aggregate  during  the  tests  in  the  three  operating 
regimes  are  given  below. 


A 

B 

C 

Flowrate 

3.17 

2.64 

2.47 

Pressure  in  the  pressurized  lines,  mH20 

4.6 

5.2 

5.4 

Electric  motor  current  load,  A 

185 

160 

150 

Maximum  amplitude  of  housing  vibration,  2 A,  //m 

60 

35 

no 

As  a  result  of  the  dynamic  studies  of  the  pump  sets  we  found  not  only  the  maximum  vibration  in  the 
elements,  but  also  found  the  frequency  range  of  vibrational  loads.  It  turned  out  that  these  vibrational 
loads  were  characteristic  ones  in  the  frequency  spectrum:  /i  =  2  Hz,  /2  =  2.5  Hz,  =  4^  Hz,  f  ^  —  5.33  Hz 
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(/rot)’  /s  —  10*66  Hz  (2/j.qJ,  and  /g  —  31.98  (/biade)»  where  is  the  rotation  frequency  and  /blade 
blade  frequency. 

Fig.  2  shows  a  graph  of  the  variation  in  peak-to-peak  amplitude  of  vibrations  (2A)  at  the  characteristic 
frequencies  as  a  function  of  the  operating  regime  of  the  pump.  We  note  that  in  regime  A,  i.e.  an  increase 
in  flowrate,  a  maximum  in  the  vibration  of  the  lower  guide  bearing  was  recorded  at  /  —  and 
f  —  f blade*  Vibration  at  lower  frequencies  were  considerably  less.  In  the  optimal  flowrate  zone  (regime 
B)  vibration  was  registered  to  be  at  the  same  level  for  almost  all  of  the  characteristic  frequencies  with 
some  increase  in  amplitude  at  /  =  2  Hz  and  f  =  2.5  Hz.  For  pump  operation  in  regime  C,  i.e.  decreasing 
delivery ,  there  is  a  redistribution  of  the  dominant  frequencies  in  the  spectrum.  The  greatest  vibration  is 
noted  at  low  frequencies  below  the  rotation  frequency.  Thus  when  the  pump  is  in  operation  with  a 
flowrate  of  less  than  Q  =  2.64  m  /s  the  possibility  of  low  frequency  vibration  (2  to  4  Hz)  goes  up  sharply, 
which  makes  the  major  contribution  to  the  integral  dynamic  load  acting  on  the  housing  and  support 
elements  of  the  pump  set,  which  may  be  a  cause  of  equipment  breakdowns  and  failures. 

Long-term  observation  of  pump  operation  has  shown  that  when  pump  delivery  is  reduced  to  that 
corresponding  to  a  current  draw  of  the  motor  of  less  than  150  A,  the  vibration,  noise  and  level  of 
hydrodynamic  load  increase  noticeably. 

We  may  note  that,  on  the  basis  of  spectral  analysis  of  the  hydrodynamic  pressure  pulsations  at  the  walls 
of  the  spiral  case,  in  all  the  studied  regimes  most  of  the  pulsational  load  energy  is  transmitted  at  those 
frequencies  for  which  an  increase  in  pump  housing  vibration  was  recorded.  In  addition,  a  significant  part 

of  the  energy  was  transmitted  at  triple  the 
rotational  frequency. 

Fig.  3  shows  the  spectral  characteristics  of 
pressure  pulsations  in  transits  of  the  spiral  case 
near  the  teeth  and  at  the  discharge  branch. 

The  power  of  pressure  pulsations  in  the  middle 
part  of  the  spiral  case  and  in  its  discharge  section 
differ  significantly  from  each  other.  These 
differences  are  especially  noticeable  in  Q  =  2,47 
m^/s  regimes  where  the  pressure  pulsation  power 
at  the  discharge  is  double  or  more  that  in  the 
middle  sections.  As  the  flowrate  increases  this 
difference  decreases.  At  flowrates  of  Q  =  2.47 
m^/s  most  of  the  pressure  pulsation  power 
increase  occurs  at  frequencies  of  2  and  2.5  Hz,  i.e. 
in  the  low-frequency  region  where  the  process  may 
be  classified  as  random.  As  the  flowrate  itself 
increases  a  rise  in  hydrodynamic  load  is  registered 
at  the  determined  frequencies  of  and  2/^.^^; 
the  maximum  contribution  occurs  at  /  =  2/^,^^. 
The  spectrum  of  pressure  pulsations  on  a 
discharge  test  bench  and  the  spectrum  of  pump 
housing  and  drive  motor  support  structure 
vibrations  are  roughly  similar;  this  indicates  that 

Figure  3.  Pressure  Pulsation  Spectra  at  the  Pump 
Discharge  in  A,  B  and  C  regimes 
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the  main  reason  for  the  increase  in  pump  vibrations  may  be  the  hydrodynamic  load  occasioned  by  the 
change  in  pump  operating  regimes,  in  which  the  regime  point  on  the  characteristic  curve  approaches  the 
boundary  of  the  performance  range.  This  in  turn  is  occasioned  by  a  change  in  the  operating  regime  of  the 
hydraulic  network  that  includes  the  given  pump. 

The  data  obtained  from  this  investigation  allows  us  to  verify  the  performance  range  of  a  pump  when 
installing  it  and  hooking  it  into  the  network,  when  the  dynamic  load  at  the  performance  boundaries  needs 
to  be  determined^.  It  is  a  good  idea  to  determine  the  coupling  between  pressure  pulsations  in  flow 
sections  of  the  pump  set  and  vibration  in  the  equipment  outside  the  performance  boundaries,  as  well  as  to 
find  those  frequencies  at  which  dangerous  increases  in  vibrations  can  arise.  It  is  clear  that  this  work 
should  be  carried  out  at  the  design  stage  of  a  PWSS  and  during  startup  testing  of  the  pumps. 

One  of  the  most  important  network  parameters  is  flowrate.  A  precise  determination  of  flowrate  for  a 
comparatively  mildly-sloping  characteristic  curve  is  one  of  the  necessary  conditions  for  finding  the  actual 
performance  boundaries  on  the  characteristic  curve.  Since  the  hard  part  in  precisely  determining  the 
location  of  the  regime  point  on  a  mildly-sloping  characteristic  curve  is  a  precise  determination  of  the 
flowrate,  and  at  present  this  task  is  fraught  with  difficulties  owing  to  the  lack  of  reliable  and  accurate 
means  of  measurement,  we  suggest  measuring  the  performance  range  from  the  amount  of  hydrodynamic 
load  acting  in  the  flow  sections  of  a  pump.  For  this  purpose,  in  the  course  of  pre-startup  testing,  control 
measurements  of  the  dynamic  loads  in  various  operating  regimes  should  be  made,  i.e.  for  various  locations 
of  the  regime  point  on  the  characteristic  curve.  For  the  most  reliable  results  when  there  are  no  reliable 
means  of  measuring  the  flowrate  for  the  operating  regime  of  a  pump  set,  the  current  draw  of  the  electric 
motor  drive  needs  to  be  monitored. 

Conclusions 

1.  In  calculating  and  designing  a  hydraulic  network  for  the  process  water  supply  system  of  an  AES  we 
need  to  account  for  the  hydrodynamic  loads  arising  in  the  flow  sections  of  pump  sets  upon  variation  in 
the  operating  regimes  of  the  consumers. 

2.  When  centrifugal  pump  sets  operate  in  a  PWSS,  the  location  of  the  regime  point  on  the  pump 
characteristic  curve  needs  to  be  monitored.  Excursions  of  the  regime  point  beyond  performance  limits 
lead  to  significant  increases  in  dynamic  load:  pressure  pulsations  and  vibration. 

3.  The  primary  reason  for  increased  vibration  of  pump  sets  is  the  effect  of  hydrodynamic  pressure 
pulsations  caused  by  the  regime  point  going  beyond  performance  limits.  The  frequency  spectra  of  pressure 
pulsations  and  vibrations  for  the  primary  assemblies  of  a  pump  set  are  for  the  most  part  identical.  As  the 
flowrate  increases  an  increase  in  the  level  of  dynamic  load  at  the  rotation  frequency  and  multiples  thereof 
is  noted.  As  the  flowrate  decreases  there  is  an  increase  in  the  dynamic  load  registered  at  low-frequencies 
in  the  range  from  2  to  4  Hz. 
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TURBINES,  ENGINES,  PROPULSION  SYSTEMS 


UDC  621.165.53 

Characteristics  of  Idle  Running  of  Turbine  Stages  and  Sections 

907F0302B  Minsk  IZVESTIYA  AKADEMII  NAUK  BSSR  SERIYA  FIZIKO-ENERGETICHESKIKH 
NAUK  in  Russian  No  1 ,  Jan-Mar  90  (manuscript  received  20  May  88)  pp  22-26 

[Article  by  G.  A.  Shapiro,  A.  D.  Kachan,  and  Ye.  I.  Efros,  Kirov  Polytechnic 
Institute  and  Belorussian  Polytechnic  Institute}  first  paragraph  is  verbatim 
IZVESTIYA  AKADEMII  NAUK  BSSR  SERIYA  FIZIKO-ENERGETICHESKIKH  NAUK  English 
abstract] 

[Text]  The  experimental  results  obtained  have  been  employed  to  derive 
generalized  relations  for  evaluation  of  the  volume  vapor  flow  rate  and  thermal 
gradient  of  turbine  stages  and  sections  with  different  geometric 
characteristics  under  zero  efficiency  conditions. 


Investigation  of  the  low— flow  modes  of  turbine  stages  is  of  great  importance 
in  increasing  the  reliability  and  efficiency  of  turbines  operating  at  partial 
loads.  Specifically,  such  modes  are  characteristic  for  the  low-pressure 
section  of  direct  heating  turbine  plants. 

Other  research^  ®  has  shown  that,  under  actual  turbine  operating  conditions, 
the  economy  of  both  the  individual  turbine  stage  and  section  with  a  specified 
geometry  virtually  singularly  determine  (given  different  combinations  of  mass 
steam  flow  and  backpressure)  either  the  magnitude  of  the  available  heat 
differential  Ho  or  the  volumetric  flow  of  steam  at  the  outlet  Gvz.  A 
reduction  in  Gvz  (or  in  Ho)  leads  to  a  reduction  in  the  relative  internal 
efficiency  oi  (and  the  used  heat  differential  Hi)  and  a  transition  of  the 
stages  from  the  range  of  active  modes  with  Hi  >  0  to  the  region  of  energy 
consumption  with  Hi  <  0.  The  idle  run,  i.e.,  that  mode  where  the  integral 
value  of  the  internal  power  of  a  stage  (section)  has  a  zero  value,  is  thus  one 
of  the  singular  operating  modes  of  a  turbine  stage  or  section.  Corresponding 
to  the  idle  run  are  specific  values  of  the  volumetric  flow  of  steam  (Gv2)x.x 
and  available  heat  differential  Hox.x  for  each  stage  (section).  Knowing  the 
values  of  (Gv2)x.x  and  Hox.x  is  of  great  interest.  This  is  connected  with  the 
fact  that  the  methodological  principles  of  determining  sections*  capacities 
changes  in  the  range  of  modes  with  positive  and  negative  efficiencies.  In  the 
range  where  a  group  of  stages  has  a  positive  power,  their  efficiency  may  be 
satisfactorily  described  by  the  relative  magnitude  of  the  reduced  ratio  of 
their  velocities,®’^  whereas  when  hoi  <  0,  the  stages'  power  requirement  is 
determined  primarily  by  the  steam  density.^ 
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Figure  1  presents  the  dependencies  of  the  relative  used  and  available  heat 
differentials  (Hi/Hih  and  Ho/Hoh),  where  the  subscript  ”H”  refers  the  the  mode 
with  the  rated  volumetric  flow  of  steam  (Gv2)h,  from  6v2  =  Gv2/(Gv2)h  for  the 
low-pressure  sections  of  T-50-130,  PT-60-130,  PT-135-130,  and  T-180-130 
turbines  during  operation  with  open  grid  valves.  The  experimental 
characteristics  have  been  obtained  on  the  basis  of  results  of  full-scale 
research,  the  method  of  which  has  been  examined  in  rather  great  detail 
elsewhere.^*®  Figures  1  and  2  also  present  data  from  the  All-Union  Heat 
Engineering  Institute  imeni  F.  E.  Dzerzhinskiy  regarding  the  last  states  of 
K-300-240  turbines  produced  by  the  Kharkov  Turbogenerator  Plant  imeni  S.  M. 
Kirov  [KhGTZ]  and  Leningrad  Metal  Plant  imeni  22nd  CPSU  Congress  [LMZ].2.3 


Figure  1 .  Dependence  of  the 
different  turbine  sections  and 
stages  used  and  the  available 
heat  differentials  on  the 
volumetric  flow  of  steam:  1-4, 
low-pressure  sections  of  PT-60- 
130,  PT-135-130,  T-50-130,  and 
T-180-130  turbines;  b,  last 
stages  of  K-300-240  LMZ  and  K- 
300-240  KhGTZ  turbines. 
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Figure  2.  Magnitudes  of  the  volumetric 
flow  and  available  heat  differential  of 
the  different  states  and  sections  in  an 
idle  run  mode:  1,  low-pressure  section 
of  a  PT-60-130  turbine;  4  through  6, 
low-pressure  section,  the  section 
consisting  of  stages  23-24,  and  stage 
23  of  a  PT-135-130  turbine;  7  and  8, 
low-pressure  section  and  stage  24  of  a 
T-50-130  turbine;  9  and  10,  low- 
pressure  section  (one  flow)  and  stage 
26  of  a  T-180-130  turbine;  and  11 
through  14,  last  stages  of  K-300-240 
LMZ,  K-300-240  KhGTZ,  and  K- 100-90 
turbines  before  and  after  updating. 


The  quantity  (Gv2)x.x  (Figu£e  1)  has  very  different  values  for  different  flow 
areas  (the  same  applies  to  Hox.x).  The  magnitude  of  the  available  heat 
differential  in  the  rated  mode  Hoh  (Figure  2)  has  a  significant  effect  on 
(Gv2)x.x;  however,  there  is  no  singular  correspondence  between  (Gv2)x.x  and 


15 


Hoh  (or  between  Hox.x  and  Hoh).  Figure  2  presents  additional  data  regarding 
the  individual  sections  and  stages  of  the  low-pressure  sections  of  the 
turbines  under  investigation.  These  data  have  been  obtained  by  using  power 
characteristics  based  on  the  method  presented  elsewhere.*  In  addition, 
extrapolation  of  the  results  of  research  conducted  at  the  All-Union  Heat 
Engineering  Institute  imeni  F.  E.  Dzerzhinskiy’  resulted  in  data  regarding  the 
idle  run  of  the  last  stage  of  a  K-100-90  turbine  before  and  after  updating. 

To  select  the  possible  criterion  possessing  the  idle  run  characteristics  of 
different  stages,  we  will,  as  a  first  approximation,  examine  a  one-dimensional 
model  of  the  steam  flow  under  the  hypothesis  that  the  flow  channels  of  the 
fixed  and  working  blades  are  completely  filled  by  the  flow.  For  the  specified 
conditions,  in  accordance  with  Euler's  equation,  the  zero  power  of  a  state 
will  occur  when  the  peripheral  components  of  the  absolute  flow  velocity  before 
and  after  the  impeller  are  in  equilibrium,  i.e.,  when 

(1) 

Cj  cos  a,  —  M  —  cos  [io, 


where  ci  and  ai  are  the  absolute  velocity  and  angle  at  which  the  flow  exits 
the  guide  apparatus,  wz  and  Bz  are  the  relative  velocity  and  angle  at  which 
the  flow  exits  the  rotating  blade  row,  and  u  is  the  peripheral  velocity  of  the 
impeller  at  its  average  diameter. 

Considering  that  ci  =  Gvi/Fc  and  wz  =  Gvz/Fi  (G  being  the  mass  flow  rate  of 
steam,  Fc  and  Fi  being  the  outlet  area  of  the  channels  of  the  nozzles  and 
moving  blades,  and  vi  and  vz  being  the  specific  volume  before  and  after  the 
rotating  blade  row),  we  obtain  the  following  from  (1): 

(Gt>2)x  .X  _ 

(2)  cosai  ,  fc 

\  ^2  /  X .  X  cos  P2  f  J, 

In  view  of  the  sharply  reduced  heat  differential  in  the  idle  run  mode,  the 
ratio  (vi/vz)x.x  included  in  (2)  is  little  different  from  unity  and  has  close 
values  for  the  different  stages.  The  same  may  be  said  regarding  the  magnitude 
of  the  ratio  cos  ai/cos  Bz .  For  example,  for  the  stages  of  the  sections  whose 
research  results  have  been  presented  above,  cos  ai/cos  Bz  is  between  1.01  and 
1.12.  The  decisive  effect  of  the  ratio  of  the  outlet  areas  Fc/Fi  on  the 
magnitude  of  the  complex  (Gvz)x.x  cos  B2/(Fiu)  thus  follows  from  (2).  It 
should  be  noted  that  the  quantity  Fiu/cos  Bz  is  the  volumetric  flow  of  steam 
behind  the  stage  in  a  mode  with  axial  output  of  steam  from  the  rotating  blade 
row. 


In  a  first  approximation,  a  stage's  available  heat  differential  may  be 
specified  as 


Ho 


2rf  \  wi  W2  / 


where  n*  is  the  maximum  relative  internal  efficiency  when  using  the  output 
velocity. 
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In  an  idle  run  mode,  allowing  for  (1),  we  obtain 


Analyzing  expressions  (3)  makes  it  possible  to  hypothesize  that,  just  like  the 
relative  volumetric  flow  (Gv2)x.x  cos  Bz/CFiu),  the  characteristic  ratio 
u2/(2Hox .X n* )  depends  primarily  on  Fc/Fi.  In  the  case  where  the  turbine 
stages  have  variable  operating  modes,  the  available  heat  differential  is 
virtually  singularly  determined  by  the  quantity  Gv2.  This  gives  additional 
basis  for  using  the  general  criterion  to  generalize  data  regarding  (Gv2)x.x 
and  Hox.x.  The  quantity  n*  specifies  the  energy  losses  when  steam  flows  in 
the  channels  of  the  cascade  and  characterizes  the  aerodynamic  perfection  of 
the  flow  area.  Proceeding  from  this,  the  efficiency  n*  is  identified  as  an 
independent  parameter  of  the  stage.  The  resultant  conclusions  apply  to 
individual  stages.  It  was  proposed  that  the  following  mean  values  be  used  for 
a  section  of  turbine  stages: 


u  = 


n 


where  n  is  the  number  of  stages  in  the  section.  The  angle  62  is  specified  as 
B2  =  arcsin(Fi/FT ) ,  where  Ft  is  the  face  area  of  the  impeller. 


Considering  what  was  said  above,  the  research  results  (Figure  3)  were 
processed  accordingly.  The  experimental  data  regarding  both  individual  stages 
and  sections  of  different  turbines  (Figure  3)  are  satisfactorily  generalized 
by  the  dependencies 


xCPsp.2  ^  f(I±. 
Fji  W  \ 


;  2u2/(2/fo..xn*)=<P  T 


The  good  agreement  of  the  data  with  respect  to  (Gv2)x.x  of  the  model  stages 
with  a  high  degree  of  fanning®  and  from  a  full-scale  experiment  should  also  be 
noted. 

The  dependencies  presented  in  Figure  3  may  be  approximated  with  sufficient 
precision  (the  broken  lines)  by  the  functions 

(4)  (0V2)x.x  PP®  p2/(^JI^)^  c/f.i) 

(5) 
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where  the  coefficients  are  as  follows:  a  =  0.64  and  b  =  2.28. 


Figure  3.  Effect  of  the  ratio  of  the  outlet 
areas  of  the  nozzle  and  rotating  blade  row  on 
the  idle  run  characteristics  of  the  following 
sections  and  stages:  1  through  3,  the  low- 
pressure  section,  sections  consisting  of  stages 
27  through  29  and  27  through  28  of  a  PT-60-130 
turbine;  4  through  6,  the  low-pressure  section, 
the  section  consisting  of  stages  23  and  24,  and 
stage  23  of  a  PT-135-130  turbine;  7  and  8,  low- 
pressure  section  and  stage  24  of  a  T-50-130;  9 
and  10,  low-pressure  section  (one  flow)  and 
stage  26  of  a  T-180-130  turbine;  11  through  14, 
last  stages  of  K-300-240  LMZ,  K-300-240  KhGTZ, 
and  K-100-90  turbines  before  and  after  updating; 
15,  model  stages’;  and  16,  an  approximation. 
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In  accordance  with  (4)  and  (5),  the  idle  run  mode  of  a  stage  or  section  occurs 
when  the  volumetric  flow  of  steam  is  reduced  to  the  value 

(Gi'2)x.x  »  auFjcos  P2: 

in  which  case  the  available  heat  differential  amounts  to 

The  relative  means  square  error  in  approximating  the  experimental  data  by 
functions  (4)  and  (5)  is  at  a  level  of  4  percent,  which  is  commensurate  with 
the  error  of  the  experiment.  Considering  this  and  bearing  in  mind  the  rather 
broad  range  of  geometric  characteristics  of  the  stages  studied,  it  is  possible 
to  recommend  dependencies  (4)  and  (5)  for  estimating  the  idle  run  parameters 
of  different  turbine  stages  and  sections. 

Conclusions 

1.  It  was  discovered  from  an  analysis  of  the  physical  processes  in  a  turbine 
stage  and  confirmed  experimentally  that^the  ratio  of  the  minimum  areas  of  the 
nozzle  channels  and  rotating  blade  row  Fc/Fi  is  a  convenient  design  parameter 
for  generalizing  the  idle  run  characteristics  of  turbine  stages  and  sections. 

2.  Generalized  dependencies  have  been  derived  that  make  it  possible  to 
determine  the  volumetric  flow  of  steam  (Gv2)x.x  and  the  available  heat 
differential  Hox.x  in  a  mode  with  zero  efficiency  for  turbine  stages  and 
sections  with  different  design  characteristics. 
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[Article  by  V.  D.  Dzhunushaliyev  and  E.  X.  Chokoyev,  Physics  Institute,  AN 
Kirghiz  SSR,  and  Kirghiz  State  University  imeni  50th  Anniversary  of  the  USSR] 

[Text]  Laser  technology  is  becoming  increasingly  important  in  modern  methods 
of  treating  materials.  This  is  precisely  why  a  great  deal  of  attention  is 
being  paid  to  basic  and  applied  research  on  this,  particularly  to  methods  of 
laser  hardening  a  metal  surface.  The  majority  of  works  in  this  direction 
relate  to  laser  treatment  in  gaseous  media,  although  the  idea  of  the 
possibility  of  laser  hardening  in  liquid  media  has  recently  been  expressed.^ 
Because  of  their  great  heat  release,  the  cooling  processes  in  these  media 
occur  much  more  intensively  than  in  a  gaseous  medium.  Treating  the  surface  of 
a  steel  product  in  a  liquid  medium  by  using  a  moving  laser  beam  increases  the 
hardness  of  its  surface  significantly. 

The  nature  of  the  occurrence  of  thermophysical  processes  in  metal  when  it  is 
treated  in  a  liquid  medium  by  a  moving  laser  beam  may  be  explained  by  solving 
the  following  heat  conduction  equation: 


(1) 


cp 


dT(x,y,z,t) 

dt 


=\M(x,y,z,t) 


here  c,  p,  and  A  are  respectively  the  heat  capacity,  density,  and  temperature 
conduction  coefficient  of  the  material,  and  T  is  its  temperature.  The 
equation  ignores  the  dependence  of  the  temperature  conduction  coefficient  on 
temperature.  Let  the  coordinate  plane  XY  lie  on  the  surface  of  the  specimen 
and  the  OZ  axis  be  perpendicular  to  this  surface.  To  simplify  the  problem,  we 
will  examine  a  rectangular  plate.  The  boundary  conditions  have  the  following 
form: 


(2) 


dr  I 

dz 


z=0 


^q^expl 


here  the  first  term  in  the  right  part  describes  the  heating  of  the  plate  by  a 
laser  beam,  qo  is  the  intensity  of  the  laser  radiation  at  the  center  of  the 
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irradiated  spot,  xo  and  yo  are  the  coordinates  of  the  initial  point  of  the 
beam's  motion,  a  is  the  characteristic  size  of  the  beam,  Vx  and  Vy  are 
projections  of  the  velocity  vector  of  the  laser  beam's  motion  in  the  plane  XY, 
and  the  second  term  describes  the  cooling  of  the  heated  plate  in  the 
surrounding  space  through  its  upper  surface. 


(3) 

dz 

(4) 

du 

Boundary  conditions  (3)  and  (4)  describe  heat  release  to  the  surrounding  space 
through  the  plate's  lower  and  side  surfaces,  respectively.  We  will  further 
assume  that  the  liquid  medium  is  water. 

It  should  be  noted  that  cooling  in  water  at  high  temperatures  is  essentially 
nonlinear. 2. 3  This  is  related  to  the  fact  that  the  following  modes  of  heat 
release  into  water  are  observed  as  a  function  of  T: 

1 .  Convective  heat  transfer 

(5)  q  ~  qconv  —  “0[conv(T  ~  To ) J  T  5  Ti 

here  qconv  is  the  coefficient  of  convective  heat  transfer,  and  To  is  the 
temperature  of  the  environment. 

2.  Nucleate  boiling 

(6)  q  =  qnucl  =  -[I3(T  -  Ts)]";  Tl  <  T  <  T2 

Ts  is  the  temperature  of  the  boiling  liquid,  and  B  and  n  are  eB^>irical 

coefficients. 2 

3.  Film-type  boiling 

(7)  q  =  qfii»  =  -afii»(T  -  Ts);  T  i  Ts 

afii»  is  the  coefficient  of  heat  release  in  the  case  of  film- type  boiling. 

4.  An  intermediate  mode  in  which  T2  <  T  <  Ts.  The  dependence  q(T)  in  this 
case  is  determined  by  the  condition  of  the  joining  of  the  modes  of  film-type 
and  nucleate  boiling. 

We  determined  the  temperatures  Ti ,  T2 ,  and  Ts  at  which  the  change  in  cooling 
modes  occurs  from  the  experimental  dependence  of  the  heat  flow  on  the 
temperature  of  the  body's  surface^  (see  Figure  1).  As  is  evident,  linear  heat 
conduction  equation  (1)  with  the  nonlinear  boundary  conditions  (2)  through  (7) 
is  practically  impossible  to  solve  analytically.  Even  solving  it  numerically 
in  full  measure  is  a  very  complicated  task.  Therefore,  to  test  the  model  in  a 
first  approximation,  we  reduced  the  three-dimensional  non-steady-state  heat 
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conduction  equation  (1)  to  a  one-dimensional  non~steady“State  equation. 
Physically,  this  means  that  the  heating  of  a  thin  rod  by  a  moving  laser  beam 
was  examined. 


Figure  1.  Piecewise-linear 
approximation  of  heat  release  by  a 
heated  surface  into  a  liquid  medium. 

In  view  of  the  fact  that  the  rod  was 
considered  infinitely  thin,  the 
processes  of  both  the  arrival  of  the 
heat  and  the  cooling  were  described 
by  including  the  heat  sources  in  the 
heat  conduction  equation: 


(8) 


dT  ,  d^T 
dt  dx^ 


-f<7v 


here  qv  is  the  right  part  of  boundary  condition  (2).  Heat  transfer  with  the 
environment  in  accordance  with  one  of  the  modes  (3)  through  (7)  was  specified 
at  the  ends  of  the  rod  depending  on  their  temperature. 


The  numerical  solution  was  performed  by  the  sweep  method  in  accordance  with  a 
scheme  described  elsewhere. Its  analysis  showed  that,  in  a  time  on  the  order 
of  several  seconds,  the  heat  transfer  process  stabilizes  and  movement  of  the 
hot  spot  (which  remains  practically  unchanged  in  shape)  behind  the  laser  beam 
is  observed  (Figure  2). 
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Figure  2.  Change  in  the 
maximum  temperature  in  the 
center  of  a  hot  spot  over 
time. 


Figure  3.  Temperature  profiles  (1,  2). 
Ratio  of  the  speed  of  cooling  in  water 
to  the  speed  of  cooling  in  air  (3);  4, 
indicates  a  laser  beam. 


Figure  3  presents  the  characteristic  temperature  distributions  that  move  at 
the  constant  rate  V.  Curves  1  and  2  describe  the  heating  and  cooling  in  air 
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and  in  water,  respectively,  given  the  following  parameters:  laser  beam  power, 
Q  =  IkW;  beam  velocity,  V  =  1  cm/s.  It  is  evident  that  under  the  effect  of  a 
moving  laser  beam,  the  rate  of  cooling  in  a  liquid  medium  is  approximately  two 
orders  of  magnitude  greater  than  the  rate  of  cooling  in  air. 

The  calculations  showed  that  in  the  rate  interval  investigated  by  us  from  0.1 
to  1.0  cm/s  and  with  a  laser  beam  diameter  of  0.5  to  2.0  mm,  the  size  of  the 
hot  spot  (defined  as  the  region  in  which  the  temperature  exceeded  half  the 
maximum  value  of  the  temperature)  depends,  for  practical  purposes,  solely  on 
the  power  of  the  laser  radiation. 

It  is  evident  from  an  analysis  of  the  solution  to  the  equation  for  heat 
conduction  at  different  speeds,  laser  radiation  powers,  and  cooling 
capabilities  on  the  part  of  the  environment  that  the  existing  asymmetry  of  the 
instantaneous  temperature  distribution  in  the  rod  is  strongly  dependent  on  the 
cooling  capability  and,  to  a  lesser  degree,  dependent  on  the  remaining 
parameters.  We  will  also  note  that  the  problem  of  the  heating  and  cooling  of 
a  specimen  hardened  by  a  moving  energy  source  has  already  been  examined. 

These  solutions  do  not,  however,  consider  cooling,  nonlinearity  in  the 
boundary  conditions,  and  limited  dimensions  of  the  heated  body. 

In  conclusion,  the  authors  express  their  gratitude  to  Zh.  Zh.  Zheyenbayev, 
corresponding  member  of  the  Kirghiz  SSR  Academy  of  Sciences,  for  his  useful 
discussions. 
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paragraphs  are  IZVESTIYA  AKADEMII  NAUK  ARMYANSKOY  SSR:  SERIYA  TEKHNICHESKIKH 
NAUK  abstract] 

[Text]  This  article  describes  the  development  and  recommendation  of  a  new 
mock-up  of  the  DVSD-1200  light  range  finder  that  includes  the  principal 
advantages  of  a  microwave  light  range  finder  under  conditions  of  strong 
perturbing  effects:  sharp  temperature  differentials,  atmospheric  turbulence, 
air  pollution,  etc.  Thanks  to  the  introduction  of  an  electronic  counter,  the 
time  required  to  make  measurements  was  reduced  by  a  factor  of  five,  which 
makes  it  possible  to  increase  the  number  of  measurements  and  reduce  the  random 
error  in  making  the  final  measurement  of  the  phase  cycle. 

Laboratory  and  industrial  tests  of  the  mock-up  on  comparators  at  the  Yerevan 
Physics  Institute,  the  Institute  of  High-Energy  Physics,  and  in  the  network  of 
the  USSR  Academy  of  Sciences'  Institute  of  Nuclear  Research  ranging  network 
showed  the  reliability  and  high-precision  of  the  mock-up's  measurement. 

Figures  3,  references  4. 

Among  research  on  the  key  parameters  and  shortcomings  of  mock-ups  of  DVSD-1200 
light  range  finders,  two  works^>2  contain  the  most  complete  information.  As  a 
result  of  the  research  conducted,  these  works  reveal  a  number  of  significant 
flaws  in  the  mock-ups:  the  limitation  of  the  temperature  range  of  the  light 
range  finder's  operation;  the  error  in  the  readout  device  making  the  final 
measurement  of  the  phase  cycle;  and  the  inadequate  stability  of  the  microwave 
oscillator's  frequency,  without  the  removal  of  which  the  production  of  DVSD- 
1200  instruments  will  not  be  feasible. ^ 

Figure  1  is  a  block  diagram  of  the  new  mock-up  of  the  DVSD-1200  light  range 
finder.  The  plane-polarized  laser  beam  (1)  of  a  type  LGN-207  A  laser  is 
directed  by  mirrors  (2)  to  a  light  modulator  (3)  on  a  KDR  [not  further 
identified]  crystal  mounted  on  a  cylindrical  cavity  with  a  central  conductor 
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1. 


(4).  The  cavity  is  mounted  on  a  positioning  table  (5)  mounted  on  guides  to 
move  the  light  modulator  140  to  150  mm  along  the  path  of  the  beam's  passage. 

A  type  Eoi  wave  is  excited  by  a  microwave  oscillator  in  the  light  modulator's 
cavity.  The  modulated  beam  is  directed  to  a  mirror-and-lens  reflector 
consisting  of  movable  objects  (7)  and  a  fixed  mirror  (7')  with  a  diameter  of 
90  mm.  A  groove  on  the  crystal  from  the  side  of  the  maximum  intensity  of  the 
cavity's  electrical  field  reduces  the  modulator's  losses  and  prevents  rays 
reflected  from  the  crystal's  outlet  face  in  the  transmission  channel  frcnn 
falling  onto  the  photoreceiver  (9)  that  is  mounted  after  the  analyzer  (8). 
Depending  on  the  position  occupied  by  the  light  modulator,  light  with  an 
intensity  of  Im  arrives  at  the  type  FEU-68  photoreceiver.  Its  amplitude 
depends  on  the  magnitude  of  the  voltage  U  on  the  crystal  and  the  intensity  of 
the  receiving  beam  I* 

(1)  4.  =  0,5/ r*  -  h.  cos i2T.DI).u)]], 

where  lo  is  a  zero-order  Bessel  function  and  Un  is  the  voltage  that  must  be 
applied  to  the  crystal  to  turn  the  polarization  90°, 

The  effectiveness  of  increasing  the  power  lies  in  the  fact  that  there  is  a 
remodulation  of  the  light,  the  maxima  of  the  intensity  of  the  light  expand  or 
double,  and  in  the  region  of  the  minima  of  the  light  there  is  an  increase  in 
the  curvature  of  the  demodulation  curves,  which  are  shown  in  Figure  2.  The 
designation  "indication  level"  in  Figure  2  refers  to  the  minimum  ratio  I»/I 
conditioning  the  differentiating  signal  on  the  light  range  finder's  indicator. 
The  experimentally  determined  "indication  level"  for  the  new  mock-up  amounted 
to  Im/I  =  1.2*10"^.  The  following  equality  takes  place  for  all  points  D'  at 
the  indication  level  according  to  which  the  position  of  the  light's  minimum  is 
determined: 

1,2-10-'  ----  0,5  { 1  -  /„  \2-UIUr.  cos  (2-D7Xm)J}. 
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V 


from  where 


To  determine  the  degree  of  the  effect  that  the  quantities  included  in 
expression  (2)  have  on  the  curvature  of  the  demodulation  curves  when  the 
indication  level  is  specified  as  I«/I  =  1.2- 10"*,  it  is  best  to  differentiate 
expression  (1)  with  respect  to  D': 


Calculations  based  on  expression  (3)  at  different  intensities  of  receiving 
light  I  and  different  ratios  U/Un  showed  that,  when  the  magnitude  of  I  changes 
by  a  factor  of  n,  the  curvature  changes  by  a  factor  of  n°-5.  In  addition,  an 
increase  in  the  light's  intensity  does  not  always  lead  to  an  increase  in  the 
measurement  precision.  The  characteristics  of  the  multiplier  phototube  show 
that,  to  increase  the  curvature  of  the  curves  in  Figure  2,  the  indication 
level  should  be  at  the  level  of  the  light's  weak  intensity.  The  pulses  from 
the  multiplier  phototube's  outlet  are  fed  to  an  amplifier  (10),  amplified  with 
respect  to  amplitude,  and  fed  to  the  signal  plates  of  a  type  6L01I  small  tube 
(11). 

Pulses  from  the  light  modulator  that  are  amplified  in  the  form  of  saw-toothed 
pulses  are  used  for  simultaneous  scanning  and  synchronization  of  the  pulses 
observed  on  the  tube's  screen,  and  the  scanning  pulses  from  the  shaper's 
outlet  (12)  are  fed  to  the  tube  (11).  In  this  case  the  scanning  amplitude  on 
the  tube  screen  shows  the  resonance  tuning  of  the  light  modulator. 

Figure  3  shows  indication  signals  that  are  characteristic  for  the  following 
operating  conditions:  a)  the  receiving  beam  has  no  fluctuations  or 
polarization  irregularities,  and  there  are  no  harmonics  in  the  modulation 
frequency;  b)  the  receiving  beam  has  a  polarization  irregularity  and 
fluctuations,  and  there  are  no  weak  amplitude  harmonics  in  the  frequency;  c) 
the  receiving  beam  is  normal,  and  there  are  weak  amplitude  harmonics  in  the 
frequency;  and  d)  the  amplitude  of  the  harmonics  is  large,  and  there  is  an 
irregularity  in  the  receiving  beam. 


26 


I 


Pulses  a  and  b  in  Figure  3  took  place  when  the  modulator  was  powered  by  a  lamp 

generator  with  an  output  power  of  P  =  100  W  without  a  isolation  gate  between 

the  modulator  and  generator.  Pulses  c  and  d  were  obtained  with  a  transistor 
generator  with  an  isolation  gate  (P  =  150  W).  The  presence  of  harmonics  in  > 

the  modulation  frequency  reduces  the  efficiency  of  light  modulation  and 

results  in  the  presence  of  residual  light. 

Two  wafers  of  KDR  crystal  were  introduced  into  the  modulator's  cavity.  One  of 
the  plates  (26)  was  connected  to  an  electric  motor  (27)  and  controlled  from 
the  instrument's  front  panel,  while  the  other  plate  (28)  was  connected  to  a 
handle  (29)  located  inside  the  instrument.  Both  plates  turn  around  their  axes 
and  change  the  modulator's  resonance  frequency.  The  handle  (29)  is  frequently 
set  in  stages  depending  on  the  temperature  of  the  medium,  while  the  electric 
motor  is  set  continuously  as  a  function  of  the  magnitude  of  the  crystal's 
heating  during  the  process  of  the  instrument's  operation. 

The  mock-up  of  the  new  light  range  finder  is  made  in  two  versions.  In  one, 
the  receiver-transmitter  has  less  of  a  slope  (12  to  14°);  in  the  other,  it  has 
more  of  a  slope  (50  to  60°).  The  first  mock-up  is  intended  for  operation  at  a 
fixed  frequency  in  a  DVSD  [not  further  identified]  mode,  whereas  the  second  is 
intended  for  operation  in  a  VSD  [not  further  identified]  mode  with  a  grid  of 
scale  frequencies.  The  basis  of  both  oscillators  is  a  GIATsINT  commercial 
thermostatted  quartz  master  oscillator  operating  at  a  frequency  of  5  MHz. 

The  set-up  for  moving  the  light  modulator  and  the  device  indicating  the  final 
measurement  of  the  phase  cycle  have  been  changed  in  the  new  mock-ups  of  the 
light  range  finder.  A  high-precision  screw  (16)  with  a  pitch  of  4  mm  that  is 
fixed  between  two  supports  is  used  for  this.  The  screw  (16)  is  connected  with 
the  light  modulator's  table  (5)  through  a  nut  (17)  made  of  fluoroplastic.  A 
friction  coupling  (19)  is  used  to  protect  the  movement  limiters  from  external 
forces  applied  to  the  handle  (18)  that  turns  the  screw  (16).  To  operate  the 
counter  (15),  a  disk  (20)  with  jaws  having  a  spacing  of  0.5  to  0.6  mm  is 
mounted  at  the  other  end  of  the  screw.  The  counting  pulses  are  shaped  when 
light  passes  from  light-emitting  diodes  (21)  through  a  slit  in  the  disk  (20) 
to  photodiodes  (22).  An  electronic  pulse  multiplier  (23)  is  used  to  shape  400 
pulses  in  one  turn  of  the  disk  (20).  The  modulator's  position  is  determined 
with  an  error  of  0.01  mm.  In  the  pulse  summer,  24  pulses  are  summed  as  the 
modulator  moves  to  one  side  and  are  removed  when  it  changes  direction.  In 
addition,  each  reading  of  the  position  of  the  minimum  of  the  receiving  signal 
is  written,  and  the  mean  value  is  output  to  the  counter's  (15)  display  screen. 
Each  reading  may  be  obtained  by  an  instruction  from  the  observer. 

The  counter  (15)  is  a  five-digit  digital  display  in  which  the  final 
measurement  of  the  phase  cycle  in  millimeters  is  illuminated,  for  example. 
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125.36.  Thanks  to  the  introduction  of  an  electronic  counter  and  readout 
summer,  the  time  required  to  make  measurements  is  reduced  by  about  a  factor  of 
5,  while  one  series  of  measurements  from  five  binary  receivers  takes  30 
seconds. 

From  a  design  standpoint,  the  new  mock-up  consists  of  three  units:  a 
receiver-transmitter,  a  microwave  oscillator,  and  a  reflector.  The  receiver- 
transmitter  is  designed  and  constructed  on  the  basis  of  separate  power  supply 
circuits  (13),  (14),  and  (25).  The  small  unit  to  supply  power  to  the  laser 
(14)  is  mounted  on  the  receiver-transmitter .  The  tube  (11),  photoreceiver 
(9),  and  amplifier  and  scanner  circuits  are  powered  from  one  source  (13), 
while  the  entire  counter  has  its  own  stabilized  power  supply  (25).  The 
receiver-transmitter's  power  requirement  amounts  to  a  total  of  17  W,  10  W  of 
which  is  used  by  the  laser.  The  transistor-based  microwave  oscillator  uses  18 
W,  8  W  of  which  are  used  by  the  quartz  oscillator's  thermostat. 

Laboratory  and  industrial  tests  of  the  mock-up  on  comparators  at  the  Yerevan 
Physics  Institute,  the  Institute  of  High-Energy  Physics,  and  the  ranging 
network  of  the  USSR  Academy  of  Sciences'  Institute  of  Nuclear  Research  in  1986 
showed  that  the  mock-up  operates  reliably,  that  its  optical  units  do  not 
require  alignment,  and  that  its  electrical  circuits  do  not  need  tuning  after 
it  has  been  transported. 
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[Text]  This  article  discusses  a  method  for  numerical  calculation  and  design 
of  the  jaw  profile  of  a  wide-range  gripping  mechanism  for  an  industrial  robot. 
The  method  makes  it  possible  to  conduct  profiling  with  any  previously 
specified  precision  depending  on  the  discreteness  selected  for  specifying  the 
component  radius. 

The  rational  significance  of  the  initial  parameter,  i.e.,  the  angle  at  the 
point  of  the  "knife  edge,"  is  examined  by  proceeding  from  the  principle  that 
imprecision  in  manufacturing  the  profile  has  a  minimum  effect  on  the  centering 
error. 

Gripping  mechanisms  are  among  the  most  important  components  of  industrial 
robots.  They  must  possess  high  flexibility  relative  to  the  parameters  of  the 
manipulation  object  and,  to  the  degree  possible,  must  not  complicate 
programming  of  the  industrial  robot.  This  is  achieved  by  the  positional 
invariance  of  the  manipulation  objects  to  their  parameters  in  the  coordinate 
system  connected  with  the  manipulator's  last  link.  For  bodies  of  revolution 
this  condition  is  met  by  using  quick-change  or  wide-range  centering  gripping 
mechanisms,  the  most  compact  designs  of  which  are  connected  with  the  use  of 
turning  jaws  with  a  special  profile. 

RTM  2  ROO-1-78  stipulates  a  method  for  the  approximate  design  of  a  jaw  profile 
along  the  arcs  of  a  circle.^  The  centering  error  (maximum  deviation  of  the 
component's  axis  from  its  nominal  position)  is  determined  by  the  dependence 

3  —  -t-  DmlnV 

(1)  “  4096A?'‘ 

where  D»ax  and  Dnin  are  the  component's  maximum  and  minimum  diameter;  R  is  the 
distance  between  the  jaw's  turn  axis  and  the  component's  axis;  and  B  is  half 
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the  angle  of  the  "knife  edge,"  i.e.,  the  angle  between  the  tangent  to  the 
component  at  the  jaw's  contact  points  (Figure  1).  It  will  be  shown  below  that 
6=01  =02.  Analysis  of  1  shows  that  the  centering  error  is  small  when  the 
changes  in  the  components'  diameters  are  small.  However,  when  Dnax  =  160  mm, 
Dmin  =0,  R  =  130  mm,  and  0  =  55°,  the  centering  error  6  amounts  to  ±0.149  mm. 
Given  the  modern  requirements  that  the  error  in  positioning  industrial  robots' 
end-effectors  be  on  the  order  of  ±0.1  mm  or  less,  such  an  error  may  be 
inadmissible.  It  is  therefore  advisable  to  develop  a  method  of  calculating 
and  constructing  a  more  exact  profile  for  centering  turning  jaws  given  an 
expanded  range  of  change  of  the  components'  diameters. 

A  method  for  numerical  calculation  and  design  of  the  profile  of  the  centering 
turning  jaws  of  wide-range  gripping  mechanisms  is  proposed. 

The  source  data  for  calculation  of  the  profile  are  as  follows:  0i ,  half  the 
angle  at  the  point  of  the  "knife  edge";  r«ax ,  the  maximum  radius  of  the 
component  manipulated;  and  R,  the  distance  between  the  jaw's  turning  axis  and 
the  component's  axis  (Figure  1).  In  accordance  with  the  theorem  of  cosines, 
from  an  analysis  of  the  triangle  ABO 

Pi=y^^+/'2— 2/?f  cosPi, 

where  Pi  is  the  radius  vector  of  the  point  B  at  which  the  upper  side  of  the 
centering  profile  makes  contact  with  the  component  and  r  is  the  component's 
radius. 

The  respective  angular  coordinate  0i  must  be 
calculated  to  fully  specify  the  position  of  the 
point  B  in  the  polar  coordinate  system  connected 
with  the  turning  jaw,  the  initial  axis  X  of  which 
passes  through  the  centering  profile  corresponding 
to  a  component  with  a  zero  radius . 

It  follows  from  Figure  1  that 

(3)  (pi  =  a— i|5i, 

where 

(4)  \l3i  =  arsin(rsinPi/pi), 
a  being  the  angular  coordinate  of  the  component 

relative  to  the  OX  axis.  Figure  1. 

When  the  increment  of  the  component's  radius  Ar  (Figure  2)  is  small,  it 
follows  from  the  triangle  A1A2D  that  A1A2  =  AaR,  the  angle  DA1A2  =  0i ,  and 
DA2  =Ar.  Consequently,  A  a  =  Ar/R  sin  0i  ,  and  when  (5)  is  integrated  with 
respect  to  r,  we  obtain 

(5)  a  =  r//?sinpi. 
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By  substituting  (4)  and  (5)  into  (3)  we  obtain  an 
expression  for  determining  0i .  In  order  to 
unequivocally  determine  the  position  of  the 
component's  axis  relative  to  the  gripping 
mechanism,  besides  the  point  B  found  on  the  upper 
side  of  the  profile,  we  must  also  fine  the 
respective  point  C  on  the  lower  side. 

Since  the  angle  02  (Figure  2)  has  not  been 
previously  specified,  it  follows  fromAAiA2E  that 
^rA2AiE  =  62.  Then,  by  analogy  with  (5), 

(6)  p2  =  acrsin(Ar/^  Aa), 

however,  in  accordance  with  (5),  Aa  =  Ar  sin  (Ji ; 
consequently. 


1,  =:  arcsin 


R  Arsin^i 
R^r 


arcsin  9, =  ^  3. 


By  analogy  with  (2)  through  (4),  it  is  possible  to  specify  the  following: 

(7)  p2=y^+''^+2i?f  C0SP2. 

(8)  11-2=  arcsin  (rsinP2/p2), 

(9)  q>2  =  a-^i2- 

For  convenience  of  use,  the  coordinates  B  and  C  may  be  transferred  to  a 
cartesian  coordinate  system: 


Xb  ~  cos  »|,  Xc  =  pj  cos  '9.,, 

(10)  ys  —  Pi  sin  9,.  yc  =  pa  sin  o.- 

Analysis  of  the  change  in  the  coordinates  of  the  points  shows  that,  as  the 
radius  of  the  component  r  increases,  there  is  initially  a  reduction  Pi,  and 
then,  from  that  moment  on,  an  increase.  Profiling  showed  that  a  segment 
constructed  when  Pi(i)  S  Pi(i-i)  is  subject  to  trimming.  It  follows  from  an 
analysis  of  (2)  for  the  extreme  that  this  phenomenon  occurs  when 

=  Rcos'^u 


Plext  =  /?Sln  t?|. 

Consequently,  rext  is  the  maximum  radius  of  the  component  at  which  a  turning 
jaw  profiled  according  to  the  method  presented  centers  the  component. 

If  the  gripping  mechanism's  design  parameters  (for  example,  the  distance  to 
that  element  of  the  gripping  mechanism  that  is  nearest  to  the  component's 


(11) 

with 

(11’) 
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axis)  allows  blank  dimensions  of  more  than  Text  and  this  capability  must  be 
used,  the  jaw  profile  may  be  adjusted  further.  All  that  is  necessary  is 
another  profiling  law  under  which  all  components  with  a  radius  greater  than 
rext  are  tangent  to  the  upper  side  of  the  profile  at  the  point  specified  for 
r  =  rext,  and  the  coordinates  of  the  points  on  the  lower  side  are  specified  in 
accordance  with  the  method  presented  below. 

It  is  assumed  that  Pi*  =  Pi  ext  =  const,  0i*  =  0iext  =  const,  xb*  =  XBext  = 
const,  and  ys*  =  ysext  =  const. 

The  angular  coordinate  of  the  point  A  of  the  component's  axis  is  determined  in 
accordance  with  the  expression 

(12)  0(0  -«(<-!)  +  Ai, 

where  a(i)  is  the  value  of  a  for  the  current  computation  cycle  and  a(i-i)  is 
the  value  for  the  previous  cycle, 

(13)  As -z 

The  value  )  is  determined  from  the  preceding  calculation  cycle,  and  that 

of  ‘Pi(i)  is  determined  from  the  expression  (Figure  1) 

^  arccos  +  P,*' -  r0/2/?p,*]. 

To  determine  2,  it  is  necessary  to  find  (on  the  basis  of  (6))  that  value  of 
B2  that  will  not  be  constant  under  the  new  conditions.  Next,  the  coordinates 
of  the  respective  point  of  the  lower  side  of  the  jaw's  profile  are  determined 
in  accordance  with  (7)  through  (9). 


By  performing  calculations  in  accordance  with  the  method  presented,  beginning 
with  r  -  0,  it  is  possible  to  determine  the  jaw  profile  with  any  previously 
specified  precision  depending  on  the  discreteness  of  the  change  in  the 
component  s  radius  that  is  selected.  Jumping  to  the  second  profiling  law  is 
possible  in  accordance  with  (11)  when  the  value  of  rext  has  been  previously 
specified.  A  conditional  jump  when  Pi  >  R  sin  Bi  may  result  in  an  error  since 
R  sin  13i  =  Pimin.  An  conditional  jump  when  Z.0BA  >  ti/2  is  more  reliable 

(15)  Pi+Tl3iCn/2. 

The  proposed  method  may  easily  be  included  in  a  simple  program  for  calculating 
the  jaw  profile  of  a  centering  gripping  mechanism  on  a  computer. 

The  angle  at  the  point  of  the  "knife  edge"  has  a  significant  effect  on  the 
shape  of  the  jaw  s  profile,  the  amount  the  jaw  opens,  the  torque  applied  to 
the  jaws,  the  dependence  of  the  centering  error  on  the  imprecision  of  the 
manufacture  of  the  jaw  profile,  and  the  reliability  of  holding  the  component. 

The  condition  of  the  reliability  of  holding  the  component  entails  having  the 
arc  distance  between  those  points  of  contact  between  the  gripping  mechanism 
and  component  that  are  most  remote  from  one  another  not  exceed  180*  -  2|i, 
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where  n  is  the  friction  angle.  For  a  symmetrical  gripping  mechanism,  this 
condition  has  the  following  form  (Figure  1):  (3i  <59°  to  69’. 

From  an  analysis  of  the  effect  of  6  on  the  torque  applied  to  the  jaws,  we 
obtain  the  following  in  a  first  approximation: 

(A:G/?/sin3)/2. 

where  K  is  the  safety  factor,  G  is  the  weight  of  the  blank,  and  f  is  the 
coefficient  of  the  component's  friction  against  the  gripping  mechanism's  jaw. 
All  these  quantities,  with  the  exception  of  B,  may  be  considered  constant. 
Consequently,  when  13  is  increased,  the  torque  increases. 

The  effect  of  the  size  of  the  angle  (3  on  the  dependence  of  the  centering  error 
on  the  imprecision  of  manufacturing  the  tooth  profile  A  is  broken  down  into 
two  cases:  the  imprecision  of  manufacture  A  on  one  side  of  the  jaw  (6i)  and 
the  imprecision  of  manufacture  A  on  both  sides  of  the  jaw  (82).  It  is  not 
difficult  to  show  that  in  the  first  case  the  centering  error  61  =  A/sin  213, 
whereas  in  the  second  case  82  =  81  +81'  or  82  =  A/sin  f3  (81  and  81'  are  the 
centering  errors  caused  by  the  imprecision  of  manufacturing  the  jaw  profile  on 
one  or  the  other  of  its  sides,  respectively). 

Analysis  of  the  constraints  on  the  angle  (3  with  respect  to  the  reliability  of 
holding  the  component  shows  that  13  5  59°  to  69*  (constraint  1-1)  (Figure  3). 
The  centering  error  in  the  presence  of  imprecision  in  the  manufacture  of  the 
profile  on  one  side  amounts  to  81  w  1.^.  To  obtain  the  same  centering 
precision  when  there  has  been  imprecision  in  manufacturing  both  sides,  it  is 
necessary  that  B  2:  45*.  The  optimum  value  of  the  angle  B  with  respect  to 
centering  precision  is  60°,  with  81  =  82  =  1.1 5A. 


t'' 


The  proposed  method  for  ntimerically  calculating  and  plotting  the  jaw  profile 
of  a  wide-range  centering  gripping  mechanism  thus  makes  it  possible  to  obtain 
a  high  centering  precision  with  a  wide  range  of  change  in  components’ 
diameters;  expand  the  ranges  of  change  of  components’  diameters,  which  is 
limited  only  by  the  gripping  mechanism’s  design  features;  and  optimize  the 


33 


gripping  mechanism's  parameters  with  respect  to  minimizing  the  effect  of  the 
imprecision  of  the  manufacture  of  the  jaw  profile,  reducing  the  force  on  the 
gripping  mechanism's  drive,  and  reducing  loose  turning  on  the  part  of  the 
jaws. 
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[Text]  This  article  examines  the  principles  lying  at  the  foundation  of  the 
formation  of  a  library  of  standardized  basic  structural  elements  for  machine 
building  components.  Algorithms  are  proposed  for  formalizing  the  processes 
examined,  and  optimization  related  to  creating  and  developing  libraries  of 
standardized  elements  is  performed. 

At  present  decisions  in  machining  process  CAD  systems  are  ii^lemented  by 
methods  of  addressation  and  synthesis.  Design  procedures  are  formulated  in 
the  process  of  examining  and  comparing  a  component  and  its  elements  with 
standardized  basic  structural  elements.  The  presence  of  a  manufacturing 
decision  that  can  only  be  used  unchanged  during  the  design  process  is  proposed 
for  each  standardized  basic  structural  element.  Design  decisions  can  only  be 
modified  and  individualized  by  using  the  synthesis  method.  Such  subsystems 
made  up  of  libraries  of  standardized  basic  structural  elements  may  be 
considered  knowledge  bases. 

The  task  of  creating  libraries  of  standardized  basic  structural  elements  and 
manufacturing  decisions  faces  all  developers  of  machining  process  CAD  systems. 
At  the  same  time  there  is  the  task  of  constantly  filling  such  libraries  with 
new  data  appearing  during  the  process  of  the  machining  process  CAD  system's 
operation  when  there  is  a  change  in  the  product  list  at  an  enterprise.  The 
correctness  of  methods  of  formulating  standardized  basic  structural  elements 
largely  determines  the  efficiency  of  the  functioning  of  the  entire  CAD  system. 
We  will  examine  the  problems  of  forming  standardized  basic  structural 
elements. 

The  set  of  components  of  an  enterprise's  product  list  may  be  represented  in 
the  form  of  two  subsets:  components  whose  production  has  already  been  set  up 
and  components  that  are  in  the  stage  of  technological  preparation  of 
production.  A  library  of  standardized  basic  structural  elements  and  a 
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corresponding  library  of 
manufacturing  decisions  is 
constructed  for  the  specified 
subsets.  Figure  1  represents  the 
process  of  the  formation  of  these 
libraries  and  shows  the 
correspondences  between  them. 

Figure  1.  Formation  of  libraries  of 
standardized  basic  structural 
elements  [YKTE]  and  manufacturing 


decisions  [TP].  (KTE  s  basic  structural  elements.) 


In  accordance  with  the  requirements  of  a  structural  systems  approach,  a 
component  may  be  described  in  terms  of  two  sets;  basic  structural  elements 
forming  a  component  and  their  structures.  In  turn,  each  basic  structural 
element  may  be  described  in  terms  of  two  subsets:  sets  of  parameters  and 
their  structural  links.  The  adequacy  of  the  model  of  a  component  constructed 
on  the  basis  of  these  subsets  is  achieved  by  selecting  the  segmentation  level 
of  the  of  a  component  into  elements  and  the  completeness  of  the  reflection  of 
information  about  these  elements  and  and  the  structural  links  between  them.^ 


To  meet  the  requirements  of  the 
structural  approach,  in  accordance 
with  a  previous  work, 2  the  basic 
structural  elements  are  divided  into 
levels  when  the  model  of  a  component 
is  constructed.  The  entire  set  of 
basic  structural  elements  forming  a 
component  may  then  be  represented  by 
joining  subsets  of  elements 
belonging  to  the  respective  levels 
of  the  hierarchy  (Figure  2). 

Figure  2.  Structure  of  the 
information  model  of  a  component 
(KTE  =  basic  structural  element; 

UKTE  =  standardized  basic  structural 


element). 


To  hasten  sorting  and  selecting  of  data  in  modern  CAD  systems,  the  coding  of 
the  basic  structural  elements  is  changed.  In  view  of  what  has  been  said,  the 
model  of  a  component  may  be  presented  in  the  form  of  the  following  expression: 

(1)  [K{E/),  KP{E/),  GKPiE/),  G(£/)|/'-=  1,...  r=:  1,  ....  /?j, 

where  Di  is  the  i-th  component  examined,  K(Etr)  is  the  code  of  the  element 
KP(Etr),  KP(Etr)  is  the  set  of  parameters  of  the  element  Etr,  G(Etr)  are  the 
structural  links  of  the  elements  Etr,  Etr  is  the  t~th  basic  structural  element 
belonging  to  the  level  r,  and  R  is  the  number  of  levels. 

In  turn,  the  basic  structural  element  of  a  component  may  be  described  in  terms 
of  the  expression 
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(2) 


E,r  ~.{K{E/),  KP{En,  QKP(E,t)\. 


For  a  standardized  basic  structural  element,  the  expression  describing  it  is 
analogous  to  expression  (2) 

(3)  {K)  =  [K{E^\  KP  OKP 

where  Eu(K)  is  the  standardized  basic  structural  element  with  the  code  K, 
KP(Eu)  is  the  set  of  parameters  of  the  element  with  the  code  K(Eu),  and 
6KP(Eu)  is  the  structure  of  the  set  of  parameters  of  the  element  Eu. 

The  set  of  all  standardized  elements  included  in  the  library  of  standardized 
basic  structural  elements  forms  the  set  of  standardized  elements  {Eu}.  This 
set  is  initially  formed  on  the  basis  of  the  set  of  components  that  have 
already  been  launched  into  production  and  is  added  to  as  new  components  are 
launched  into  production.  A  generalized  algorithm  for  forming  a  set  of 
standardized  basic  structural  elements  is  presented  below. 

Step  1.  For  elements  having  one  code,  the  subsets  that  are  common  regions  for 
both  the  parameter  sets  and  the  structure  of  the  parameter  sets  of  the 
elements  under  examination  are  found 

(4)  \KP{E^)\  rr  [KP(E,.{)\  n  -  n  \KP{E,.i)\  n  -  n  \kp{e,.a)\, 

(5)  \GKP{E^)\  =z  {G/vP(£'*.i)l  n  -  n  \GKP{E„,i)\  fl  -  fl  \GKP{Eu.a)], 

where  {l^(Ek)}  and  {6KP(Ek)}  are  sets  of  common  elements  of  the  parameter  sets 
and  the  structure  of  the  parameter  sets  formed  for  elements  (with  the  code  K) 
of  the  components  Di ,  Dz ,  ...,  Da;  {KP(Ek,i)}  is  that  set  of  parameters  with 
the  code  K  that  belong  to  the  component  Di ;  and  {GKP(Ek,i)}  is  the  set  of 
structural  links  of  the  set  of  parameters  for  those  elements  with  the  code  K 
belonging  to  the  component  Di . 

Step  2.  The  original  components  of  the  set  of  parameters  and  its  structure 
for  each  basic  structural  element  having  the  given  code  are  determined. 

(6)  \KP{E,.,)  =  {\1^{E,)\U\KP{E,4)\\[KP{E^M. 

(7)  ^GKP(Ek,^)  =  (\GKP{E„)\U\GKP[Ek,i)\\\GKP{Ek::)\, 

where  AKP(Ek,i)  and  AGKP(Ek, i)  are  original  components  of  the  set  and 
structure  of  the  parameters  for  those  elements  of  the  code  K  belonging  to  the 
component  Pi • 

Step  3.  The  previously  formed  sets  {KP(Ek,i)}  and  {GKP(Ek.i)}  of  the 
parameter  sets  for  the  standardized  basic  structural  elements  with  the  code 
under  examination  are  supplemented  with  the  resultant  data^  i*e.,  AKP(Ek,i) 
and  A6KP ( Ek , i ) 


,  _  \\KP{E^)\  U  \\KP{Fk,\)U \  KP{Ek.i)U...UKP{Ek.A)\ 

(8)  U  \^GKP{Ek,\)U  ^GI\P(Ek.2)U  ...UGKP[Ek.A)\- 
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step  4.  We  obtain  the  standardized  basic  structural  elements  with  the  given 
code  by  processing  information  regarding  the  list  of  components  manufactured 
within  the  range  of  one  code. 


Step  5.  As  a  result  of  the  formation  of  standardized  basic  structural 
elements  with  respect  to  all  codes  encountered  in  the  list  of  components,  it 
is  possible  to  produce  a  library  of  standardized  basic  structural  elements 


(9) 


'F  I  — 


[Eu,,\  U  U[Ea.x\ 


where  N  is  the  niunber  of  types  of  component  elements  throughout  the  entire 
product  list. 

When  the  product  list  changes,  a  search  is  conducted  for  new  components  of  the 
parameter  set  and  its  structure  for  each  element  code.  The  resultant  data  are 
used  to  supplement  the  information  files  of  standardized  basic  structural 
elements  describing  the  set  of  parameters  and  its  structure.  In  the  case 
where  a  new  code  is  encountered,  it  is  processed  in  accordance  with  the 
algorithm  for  forming  a  standardized  element  that  was  examined  above. 


The  algorithms  presented  make  it  possible  to  formalize  the  process  of  forming 
standardized  basic  structural  elements.  Examining  the  formalized  technique 
for  forming  them  makes  it  possible  to  optimize  the  job  of  creating  and 
developing  libraries  of  standardized  elements. 


Conclusions 


1.  A  method  of  creating  a  library  of  standardized  basic  structural  elements 
has  been  developed  on  the  basis  of  analyzing  a  component  and  the  links  between 
the  elements  forming  it. 

2.  The  distinction  of  this  method  is  that  it  affords  the  capability  of  self¬ 
organization  of  libraries  of  standardized  basic  structural  element  during  the 
process  of  their  use. 
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Computer  Calculation  of  Mode  for  Arc  Welding  Butt  Joints  With  Beveling  of 
Edges 
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[Article  by  Ye.  B.  Maslova,  student,  A.  M.  Rybachuk,  docent,  and  A.  V. 
Konovalov,  engineer;  first  paragraph  is  IZVESTIYA  VYSSHIKH  UCHEBNYKH 
ZAVEDENIY:  MASHINOSTROYENIYE  abstract] 

[Text]  This  article  proposes  an  algorithm  for  selecting  the  optimum  method  of 
welding  and  calculating  a  welding  mode  with  the  necessary  surfacing  area  and 
maximum  penetration  capability. 

When  welding  butt  joints  20  to  40  mm  thick  it  becomes  necessary  to  select  that 
welding  method  and  mode  ensuring  filling  of  the  bevel  and  complete  fusion  of 
the  cross  section.  Automatic  submerged-arc  welding  and  CO2  shielded-arc 
welding  are  most  suitable  for  products  made  of  low-carbon  steel.  Experimental 
analysis  requires  great  expenditures  of  labor,  time,  and  materials.  It  is 
proposed  that  the  task  be  performed  by  using  a  computer  and  that  subsequent 
experimentation  only  refine  the  results  obtained. 

The  work  was  performed  for  a  specific  product  —  a  cylindrical  shell  made  of 
low-carbon  steel  and  with  an  outer  diameter  of  515  mm  and  a  thickness  of  34 
mm.  The  shell  was  welded  with  a  quality  lengthwise  seam.  The  term  "quality" 
is  understood  to  mean  complete  fusion  of  the  cross  section  and  achievement  of 
the  required  fusion  area.  Owing  to  the  uneven  deformation  of  the  metal  during 
roll  forming,  the  spread  of  the  edges  in  the  joint  has  a  variable  value  for 
different  shells,  which  in  turn  makes  it  necessary  to  weld  different  amounts 
of  metal  into  the  bevel.  The  bevel  has  a  nonstandard  V-shaped  form  with  an 
angle  that  changes  along  the  length  of  the  joint  and  that  differs  in  different 
products.  The  purpose  of  this  work  is  to  select  that  welding  mode  ensuring 
filling  of  the  bevel  and  complete  fusion  of  the  cross  section  in  submerged-arc 
welding  and  CO2  shielded-arc  welding. 

To  obtain  the  source  data,  the  spread  of  the  edges  y  (Figure  1)  was  measured 
in  three  different  cross  sections  along  the  length  of  the  joint  in  35 
products.  Values  of  the  mean  square  deviation  of  the  spread  of  the  edges  of 
01  =  1.25,  02  =  1.26,  and  03  =  1.03  were  obtained  for  each  cross  section 
separately.  These  deviations  were  obtained  in  accordance  with  the  following 
formula 
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where  n  is  the  number  of  measurements  and  y  and  yr  are  the  mean  and  running 
values  of  the  spread  of  the  edges. 


Figure  1.  Joint  of  a  shell  (a)  and 
diagram  of  the  arrangement  of  the  cross 
sections  in  which  the  overlapping  of 
the  edges  were  measured. 

Q) 

For  each  cross  section,  we  find  the 

variance  in  accordance  with  the  formula 

0=0^:  ai2  =  1.57,  az^  =  1.61,  and  =  1.07. 


Estimating  the  uniformity  of  the  variances  in  accordance  with  Fisher's 
criterion  showed  that  all  of  the  values  of  the  spread  of  the  edge  for  the 
three  different  cross  sections  may  be  combined  into  one  general  set,  i.e.,  the 
variances  are  uniform.  A  mathematical  expectation  of  the  value  of  the  spread 
of  the  edges  in  the  joint  of  M(y)  =  14  mm  was  calculated  for  the  general  set. 
Statistical  processing^  showed  that  the  distribution  of  the  values  of  the 
spread  of  the  edges  is  subjected  to  a  normal  gaussian  distribution  law  with 
the  following  parameters:  a  =  14  and  o  =  1.26. 

Figure  2  presents  a  histogram  whose 
ordinate  represents  the  relative 
frequencies  h  of  the  appearance  of 
one  value  for  the  spread  of  the  edges 
or  another  in  the  grouping  intervals 
and  whose  abscissa  represents  the 
bounds  of  these  intervals.  Figure  2 
shows  a  normal  distribution  curve 
approximating  the  histogram.  The 
same  statistical  calculation  results 
may  be  obtained  by  using  a  computer. 

Figure  2.  Results  of  statistical 
processing:  a,  curve  with  a  normal 
distribution;  b,  histogram. 

A  program  to  calculate  modes  of  9 

submerged-arc  welding  and  COz 

shielded-arc  welding  in  accordance  with  the  existing  methods^  was  compiled  to 
select  the  welding  method  and  determine  that  welding  mode  providing  the 
required  depth  of  fusion.  The  constant  thermophysical  coefficients  are 
selected  from  a  previously  formulated  data  file.  Input  from  a  display  screen 
are  the  type  of  welding  (COz  shielded-arc  welding  or  submerged-arc  welding); 
the  thickness  of  the  metal  being  welded;  the  bounds  of  the  measurements  of  the 
welding  current  I,  the  welding  speed  Vb ,  the  electrode  diameter  ds,  the 
electrode  tip  b,  and  the  bevel  area;  the  surfacing  error;  and  the  measurement 
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intervals  of  the  current,  electrode  diameter,  throat  capacity,  and  welding 
speed. 

The  calculation  is  reduced  to  determining  that  welding  mode  during  which  we 
achieve  a  balance  between  the  volumes  of  the  bevel  and  surfaced  metal  given 
different  combinations  of  welding  speed,  arc  voltage,  and  electrode  diameter 
and  tip. 


Figure  3.  Flowchart  of  the  program. 
Key: 


1.  Beginning 

2.  Input  source  data 

3.  Modification  of  values  of  mode 

parameters 

4.  Determination  of  surfacing  area 

5.  Bevel  is  filled 

6.  Determination  of  depth  of  fusion  H 

7.  H  >  Hnax 

8.  Preservation  of  running  parameters 

9.  H.ax  =  H 

10.  Print  parameters  of  mode  selected 

11.  End 

12.  No 

13.  Yes 

The  calculation  algorithm  is  presented 
in  the  form  of  a  flowchart  (Figure  3). 
The  calculation  begins  with  a 
calculation  of  the  fusion  coefficient: 


(1)  ap  =  a'p  +  Av 

In  the  case  of  welding  using  a  direct  current  with  reverse  polarity 

a'p=  11,6  4  0,4  r/A-M. 

Aap  =  3600 

The  value  of  Qicit  is  calculated  in  accordance  with  the  formula 

(2)  Q„„  0.1884y>,rf2 


The  coefficients  13,  pi ,  and  p2  are  calculated  in  accordance  with  the  formulas 
(3)  0,576/2  po/cry, 


V. 


a  aQ 


_  1)  i 


« (r„,  —  7„) 

Pt 


]• 
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where  Ve  is  the  relative  speed  at  which  the  wire  is  fed  (m/hr), 

(5)  K. ---  a'py/SeOOT; 
and  j  is  the  current  density  (A/mm^) 

The  next  step  is  that  of  determining  the  surfacing  area  Ph.  We  will  estimate 
the  losses  due  to  melting  and  spattering.  For  welding  in  CO2, 

(6)  'I'„  —  4,72  -h  0,176/  —  0,000448/, 

for  submerged-arc  welding  iPp  =  0.03  (%). 

Next,  we  fine  the  rate  at  which  the  electrode  wire  is  fed  (m/hr): 

(7)  V,in  “  “p//7 

(8)  Fff~  (FyX^nnl^cn) 

where  oh  =  ap(1  -  (W*/100))  for  welding  in  CO2  and  oh  =  ap(1  -  0.03)  for 
submerged-arc  welding. 

The  surfacing  area  is  calculated  within  the  cycles  by  a  one-by-one  trial  of 
the  values  of  the  welding  current,  welding  speed,  and  diameter  and  tip  of  the 
electrode  (Figure  3).  As  soon  as  the  difference  between  the  values  of  the 
surfacing  area  and  the  bevel  area  turns  out  to  be  less  than  the  allowable 
surfacing  error  c,  the  calculation  of  the  surfacing  area  is  halted.  It  is 
advisable  to  determine  the  modes  for  the  maximum  bevel  area  value.  The 
calculation  results  in  that  range  of  combinations  of  mode  parameters  providing 
the  necessary  surfacing  area.  This  is  because  complete  depth  of  fusion  in  so 
narrow  a  bevel  is  not  possible  during  any  welding  mode.  The  fusion  capacity 
of  the  modes  was  compared  for  the  case  of  surfacing.  The  depth  of  fusion  is 
determined  in  accordance  with  the  formulas 

(9)  for  submerged-arc  welding  =  0, 1 56 

(10)  for  CO2  shielded-arc  welding  7/np  —  0,165  ]/^/V^Vnp 
where  q  =  lUnni  and  U^itp  =  k'(19  -  0.01I)(dEUD/I) . 

With  a  current  density  of  j  <  120  A/mm^ ,  k'  =  0.282j®-i’25  during  welding  with 
a  direct  current  having  direct  polarity,  and  k'  =  0.367j®i’25  during  welding 
with  a  direct  current  having  reverse  polarity.  When  j  >  120  A/mm^,  k'  =  0.92 
with  a  current  having  direct  polarity,  and  k'  =  1.12  with  a  current  having 
reverse  polarity. 

After  having  determined  the  depth  of  fusion  (in  accordance  with  (9)  and  (10)) 
for  the  possible  welding  modes,  we  select  the  one  mode  with  the  maximum  depth 
of  fusion  (Figure  3).  The  program  developed  thus  makes  it  possible  to  select 
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the  method  and  calculate  the  mode  of  welding  providing  the  required  surfacing 
area  and  the  maximum  fusion  capacity. 

Bibliography 

1.  Akulov,  A.  I.,  Belchuk,  G.  A.,  and  Demintsevich,  V.  P. ,  "Tekhnologiya  i 
oborudovaniye  svarki  plavleniyem.  Uchebnik  dlya  studentov  vtuzov"  [Technology 
and  Equipment  for  Fusion  Welding.  Textbook  for  Students  at  Higher  Technical 
Schools],  Moscow,  Mashinostroyeniye,  1977,  432  pages. 

2.  Maslov,  B.  G. ,  "Nerazrushayushchiye  methody  kontrolya  kachestva  svarnykh 
shvov"  [Nondestructive  Methods  of  Testing  Quality  of  Welds],  Moscow, 
Mashinostroyeniye,  1976,  72  pages. 

COPYRIGHT:  "Izvestiya  VUZov.  Mashinostroyeniye",  1990. 


43 


INDUSTRIAL  TECHNOLOGY,  PLANNING,  PRODUCTIVITY 


UDC  621.9 

Constructing  Matrix  Models  of  Precision  of  Manufacturing  Processes 
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[Article  by  L.  K.  Sizenov,  doctor  of  technical  sciences,  and  A.  A.  Gusev, 
candidate  of  technical  sciences;  first  paragraph  is  IZVESTIYA  VYSSHIKH 
UCHEBNYKH  ZAVEDENIY:  MASHINOSTROYENIYE  abstract] 

[Text]  This  article  examines  methods  of  constructing  multidimensional 
Diultif actoral  mathematical  models  of  the  precision  of  a  set  of  interconnected 
operations  entailed  in  manufacturing  processes  as  the  foundation  of 
forecasting  and  controlling  production  quality. 

To  model  manufacturing  processes  we  will  introduce  the  following  spaces  into 
our  examination.  The  factor  space  of  a  process  consisting  of  n  operations  is 
the  linear  space  of  the  vectors  whose  coordinates  are  the  errors  in  the 
initial  factors  of  the  blanks  and  the  parameters  of  the  manufacturing  systems 
in  which  the  operations  are  implemented.  In  accordance  with  the  linear 
dependence  of  the  factors,  it  may  be  represented  in  the  form  of  the  direct  sum 
oi  n+1  spaces  Uo  ®  Vi  ®  ...  ®Vn,  where  Uo  is  the  space  of  the  arbitrary 
vector  of  the  input  variables  of  the  initial  blanks  (error  in  dimensions, 
deviations  in  shape,  microgeometry,  hardness,  etc.);  Vi  is  the  space  of  the 
arbitrary  vector  of  the  manufacturing  system's  input  variables  (rigidity,  tool 
wear,  temperature  deformation,  adjustment  error,  etc.)  affecting  the  output 
variables  of  the  i-th  manufacturing  operation  (i  =  1,  ...,  n)  characterizing 
the  parameters  of  the  quality  of  the  blanks  undergoing  machining  (error  in 
dimensions,  deviations  of  the  shape,  arrangement  of  the  surfaces,  waviness, 
etc. ) . 

Let  the  parameters  of  the  component's  quality  that  changed  or  that  arose  at 
the  outlet  from  the  i-th  operation  form  the  arbitrary  vector  varying  in  the 
space  Ui .  This  i-th  operation  may  be  represented  as  a  linear  heterogeneous 
transform  of  the  spaces 

U;.,  ©  o  -  0  v„  -V  0  0  ...  0 

that  any  vector 

^  L  (B  ^ 
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transforms  into  the  vector 


e  Ui  o  Ki+i  e ...  0  v„ 

in  accordance  with  the  formula 

(1)  =  (^'=1 . 

where  Aio  is  the  matrix  column  of  free  terms  and  Ai  is  the  matrix  of  the 
transfer  coefficients.''  For  convenience  of  writing,  we  will  complement  the 
vector  xi  with  a  fictitious  variable  to  the  vector 


(2) 


(i  -  0,  1 . n) 


and  will  introduce  the  new  matrix 


(3) 


0  0  . 


(i 


By  taking  (2)  and  (3)  into  account,  we  may  write  transform  (1)  in  homogeneous 
form 

(4) 

Let  us  write  the  leveling  of  the  constraint  on  the  output  quality  parameters 
of  components  from  the  i-th  and  j-th  operations  ( j  <  i)  in  homogeneous  form 

(5)  =  «■  >,  =  «■' . 

or  in  heterogeneous  form 


Xi  —  +  ^{.^<-1.0  +  AiAi.iAi^2.o  +  •••  A-Ai-  ...  -Aj+tXj. 

A  manufacturing  process  may  be  represented  in  terms  of  a  graph  (Figure  1) 
illustrating  the  model  in  scalar  form  (Figure  1a)  or  in  terms  of  a 
heterogeneous  matrix  (Figure  1b  and  (6))  or  homogeneous  matrix  (Figure  1c  and 
(5)).  The  nodes  of  the  graph  shown  in  Figure  1c  that  are  indexed  by  vectors 
are  the  output  quality  parameters  of  the  process,  ^hereas  the  ribs  are 
transfer  matrices.  The  path  from  node  to  node  xi  (i  >  j)  drawn  along  the 
graph  in  the  direction  of  the  orientation  designates  the  transition  from  the 
properties  to  the  quality  parameters  iti  in  accordance  with  (5). 

By  applying  the  operation  of  mathematical  expectation  to  matrix  equation  (5), 
we  obtain 
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Figure  1.  Graph  of  a  manufacturing  process:  a,  for  scalar  writing  of  the 
constraint  equation;  b,  for  matrix  writing  of  heterogeneous  constraint 
equations  (6);  c,  for  matrix  writing  of  homogeneous  constraint  equations  (5). 


The  correlation  matrix  of  the  arbitrary  vector  of  the  output  quality 
parameters  xi  derived  as  a  result  of  linear  transformation  (5)  of  the  vector 
xj  of  the  original  manufacturing  factors  is  determined  in  accordance  with  the 
formula 

(8)  KTJ  ~  •Aj+fKTjTjAj+^-  ... ‘A/  (4  >  y  =  0, 1, ....  rt). 

For  heterogeneous  constraint  equations  (6),  we  may  write  the  following  instead 
of  the  probability  characteristics  (7)  and  (8); 

(9)  wjT  Ai^  -f-  +  ...  ..■•Aj^..Ajj^\,Q'{-A^-  ...  ■  A^+,mT. 


•  ■  Aj^.KTj  .  ...  -A/  (/  >  y  --  0.  1 ,  ....  n). 


By  isolating  the  diagonal  elements  in  correlation  matrix  (10),  we  obtain  the 
variance  of  the  output  quality  parameters  of  the  i-th  operation 
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(11) 


A/,-.-  -  ■Aj,,/<T.TjA'j.r  - 
(/r^I,2 . n,\  k  :r.(\.2,  ...,s^). 


where  Dxk<^>  is  the  variance  of  the  k-th  output  of  the  i-th  operation  and  Aik 
is  the  k-th  line  of  the  matrix  Ai. 

By  switching  in  (11)  from  the  variances  to  the  stray  fields  of  the  errors,  we 
derive  a  constraint  equation  for  the  stray  fields  of  the  output  quality 
parameters  of  the  i-th  and  j-th  operations 

(i)  •  ...  •  Aj+ikT^^Tj  X, 

^  X  Xj  X.  ^XjkXjA^ ... 

(/  >  y,  y  =r  0,  1 . rt  —  I ;  /  =;  I ,  ....  rt;  k  —  1 . iy). 

whereAxk<^>  and  kxk<‘)  are  the  stray  field  and  coefficient  of  the  respective 
scatter  of  the  k-th  output  quality  parameter  of  the  i-th  operation  andAxj  and 
kxj  are  diagonal  matrices  of  the  stray  fields  and  the  coefficients  of  the 
relative  scatter  of  the  components  of  the  vector  . 

Table  1 .  Construction  of  a  mathematical  model  of  the  precision  of  an 
automatic  line. 


(1) 

Sianw  Tex- 
Ho.iorHMec- 

'foMHOCTHbie 

xapaKTcpH- 

CTHKH 

^^^([)(J)H- 

UHCHT 

Koppe- 

i  vlHUHH 

^^^riepeiiaTOHHbie 
XapaKTepHCTHKH  [1] 

iilK^nepCHH  COOCTBCKHblX 
norpemHOCiefi 
oOpauoTKif  [1] 

Koro  npo- 
uecca 

MM 

’-Vi . 

MKM 

(M  m) 

'■-Vi  -V/-, 

’-V/-, 

■’■’-\<'  “'''Vi-Vi_i)* 
MKM- 

(  M  m^) 

SaroTOBKa 
(1=0)  (5) 
ToKapnaH 
aBTOMaTHaa 
o6pa6oTKa 

(7) 

14 

10,3 

100 

10 

0,5 

(i,u=IO,3— 0.005-14--9,G 

1 

1 

=-'^,^=lU-(l-0.5-)  =  75 

BecucHTpo- 

BO'IUJIHC^O- 

BaJibHaa 

o6pa6oTKa 

U=2) 

(8) 

9,983 

4 

0,4 

1*2=0.4-1=0,16 

rt2o=9,98.3-0, 16-10, ••)■- 
=8.3.35 

=A.v,=4-(1-0,4O=!-3,-14 

1 .  Steps  in  the  manufacturing 

process 

2.  Precision  characteristics 

3.  Correlation  coefficient 

4.  Transfer  characteristics^ 


5.  Variance  of  the  inherent 

machining  errors^ 

6.  Blank  (i  =  0) 

7.  Automatic  turning  (i  =  1) 

8.  Eccentric  grinding  (i  =  2) 
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As  an  example,  let  us  examine  the  construction  of  a  mathematical  model  of  the 
precision  of  machining  rollers  on  interconnected  corridors  of  an  automatic 
line.  Investigation  of  a  batch  of  blanks  throughout  all  stages  of  the 
manufacturing  process  resulted  in  the  precision  characteristics  mxi  and  oxi 
and  establishment  of  linear  correlation  links  rxixi-i  between  them.  Table  1 
presents  the  sequence  of  the  calculations. 

The  mathematical  model  of  the  machining  precision  is  as  follows: 


3"  -  — 0^  -f-d'-h  „  ^  *•  —  l5,36-r6,4’  I0”'’3"i*  • 

Xj  x-iXi  -vi.ro  '  -  -  JTo  ’  .\o 


Analyzing  mathematical  models  permits  well-founded  designation  of  the 
requirements  regarding  a  blank's  precision  and,  in  some  cases,  expands  the 
tolerances  for  its  manufacture.  This  in  turn  improves  the  technical  and 
economic  indicators  of  the  manufacturing  process. 


A  method  of  constructing  matrix  models  that  permits  a  comprehensive  analysis 
of  the  precision  of  an  entire  manufacturing  process  with  an  allowance  for  the 
effect  of  an  entire  set  of  factor  and  examination  of  all  of  the  principal 
quality  characteristics  of  components  has  thus  been  developed. 
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